BETHE ANSATZ AND INVERSE SCATTERING TRANSFORM 
IN A PERIODIC BOX-BALL SYSTEM 
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Abstract. We formulate the inverse scattering method for a periodic box- 
ball system and solve the initial value problem. It is done by a synthesis 
of the combinatorial Bethe ansatze at q = 1 and q = 0, which provides the 
ultradiscrete analogue of quasi-periodic solutions in soliton equations, e.g., 
action-angle variables, Jacobi varieties, period matrices and so forth. As an 
application we establish explicit formulas counting the states characterized 
by conserved quantities and the generic and fundamental period under the 
commuting family of time evolutions. 



1. Introduction 

In [1, 2], a class of periodic soliton cellular automata on one dimensional lattice 
is introduced. They are associated with crystal basis of non-exceptional quantum 
affinc algebras Uq{2„) [3]. In this paper we focus on the simplest case 0„ — ylj"'^-' 
known as the periodic box-ball system [4] , and solve the initial value problem by the 
inverse scattering method. As a result, all the properties and formulas conjectured 
in [1, 2] are established in this case. Our approach is a synthesis of the two versions 
of the combinatorial Bethe ansatz at g = 1 [5, 6] and q — [7]. We provide a self- 
contained proof for all the important statements together with several examples. 
Here is a typical time evolution pattern on the length 13 lattice: 
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Regarding the letter 1 as the background, one observes the three solitons traveling 
periodically to the right with the velocity equal to the amplitudes 3, 2 and 1 as long 
as they stay away from each other. Under the collisions they do not smash into 
pieces nor glue together. This is due to the underlying quantum group symmetry. 
The time evolution from t = 3 to t = 4 states for example have been determined 



from the following diagram: 

1122211112212 
122 -|- 1 12 -|- 1 1 1 -|- 112 -|- 122 -|- 222 -|- 122 -|- 1 12 -|- 1 1 1 -|- 111 -|- 1 12 -|- 122 -|- 112 -|- 122 
2211122211121 

Here the local vertices stand for the Uq (A^^' ) quantum R matrix sending the spin 
^ states on the NW edges to the spin states on the SE edges. This is a standard 
diagram for the row transfer matrix in the 3x1 fusion vertex model. The pecuhar 
feature here is the specialization to g = 0, which makes the R a deterministic map 
called the combinatorial R [8, 9]. The above ones correspond to the I = 3 case 
in Figure 2.1, and the associated time evolution is named T3. On the horizontal 
edge, the leftmost and the rightmost ones are the same element 122 reflecting the 
periodic boundary condition. Denote the array on the top and the bottom lines by 
p and p', respectively. It turns out that p' is uniquely determined from p by the 
combinatorial R and the boundary condition said above. The resulting map defines 
a time evolution p' = T3{p). The periodic box-ball system is a dynamical system 
endowed with such time evolutions Ti, T2, . . ., which are essentially the commuting 
family of fusion transfer matrices at g = 0. In other words, it is a solvable vertex 
model [10] at q = on the periodic lattice. It is the periodic extension of the 
original box-ball system on the infinite lattice [11]. 

In this paper we solve the initial value problem of the periodic box-ball system. 
It is done by a synthesis of the Bethe ansatze [12] at q = 1 [5, 6] and q = [7], 
which yields a periodic ultradiscretization of the inverse scattering method in soliton 
theory [13, 14]. As an application, we are able to answer the interesting question; 
what is the fundamental period of p under the time evolutions Ti,T2, . . .7 Namely, 
the minimum positive integer Ni satisfying Tf^' (p) = p. For instance for the above 
p, one has A^i = 13, A^2 = 91, = 273^. An explicit formula for the fundamental 
period under has been discovered in [15] by a combinatorial argument. Here we 
take the inverse scattering approach and establish the results (4.22) and (4.26) for 
general T;. Under our scheme, the period Af acquires the intrinsic characterization 
Afh e r (4.31) in terms of straight motions on the set 

This is an 'ultradiscrete Jacobi variety'. See (4.30). The conjectural formula for the 
generic period in the most general An"^ case [2] also admits a similar interpretation. 

Let us explain the content of the paper in some detail. In Section 2, we begin 
with the rudiments of the crystal basis theory, formulate the periodic box-ball 
system and explain its basic properties. The fundamental role is played by the 
combinatorial R. It governs the local dynamics of the system and ultimately its 
whole aspect together with the global constraint imposed by the periodic boundary 
condition. The commuting family of time evolutions Ti and the conserved quantity 
El called energy are constructed. It is done by a simple extension of the arguments 
in [16, 17, 18] to periodic versions and supplements the original description in [4]. 
Th(^ both Ti and the Ei enjoy the symmetry under the extended affine Weyl group 



It is a general feature of the system that T;(p) = Too(p) hence A'; = Noo for / > ^o, where 
lo is the amphtude of the largest soliton involved in p. See Proposition 2.5 and the subsequent 
remark. 
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') (Proposition 2.3). The time evolution Ti with sufficiently large I admits 
a simple description as a translation in (Proposition 2.5). Although quite 

straightforward, the proof of these properties announced in [1, 2] are presented here 
for the first time. 

Section 3 is a main part of the paper where we present our inverse scattering 
formalism. Wc invoke the Kerov-Kirillov-Rcshctikhin (KKR) bijcction between the 
rigged configurations and highest paths. We call it the KKR theory or combinatorial 
Bethe ansatz at g' = 1 in this paper. The necessary facts are available in Appendix 
A as well as in the original papers [5, 6, 28] and the latest review [19]. 

Roughly speaking, our action-angle variables are the rigged configurations. How- 
ever, the KKR theory works only for highest paths whereas the periodic boundary 
condition brings all the paths into the game. To reconcile them in a remarkable 
harmony via the special prescription (3.12) is the heart of our inverse scattering 
formalism. Proposition 3.7 is the crux to guarantee the well-definedness of the 
resulting direct and inverse scattering map <I> (3.15). Our solution of the initial 
value problem is presented in Theorem 3.12. It is also applicable, with a drastic 
simplification, to the original box-ball system [11] by taking the system size infinite. 

In Section 4 wc explain the origin of our inverse scattering formalism in the 
other framework of the combinatorial Bethe ansatz at g = [7] . The central role 
is played by the linear congruence equation called the string center equation (4.3). 
It is the Bethe equation on the string centers at g = whose off-diagonal solutions 
yield the weight multiplicities of the module (C'^)'^^. By the map \I' (4.18), we 
link the logarithmic branch in the string center equation with the angle variable 
in Section 3. Then our key construction (3.12) is nothing but the equivalence 
relation among the off-diagonal solutions to the string center equation (Theorem 
4.3). Through our direct and inverse scattering transforms, the nonlinear dynamics 
in the periodic box-ball system becomes a straight motion of the Bethe roots or 
angle variables. With these features (Corollary 4.4) in hand, it is straightforward 
to derive explicit formulas for the generic (4.22) and the fundamental period (4.26) 
under any time evolution Ti. The number of states characterized by conserved 
quantities (4.21) are obtained and the Bethe eigenvalue is shown to be a root of 
unity related to the generic period (Proposition 4.11). The number of disjoint orbits 
under the commuting family of time evolutions is also determined in (4.35). Our 
angle variables live in the set (4.30), which serves as an ultradiscrete analog of the 
Jacobi variety in the classical theory of quasi-periodic solutions to soliton equations 
[20, 21]. 

The results in Sections 3 and 4 uncover a significant interplay between the KKR 
theory and the string center equation. They substantially achieve a synthesis of 
the combinatorial Bethe ansatze at g = 1 and q = 0, which was already foreseen 
in the earlier works [1, 2] on generalized periodic box-ball systems. Section 5 is a 
summary. 

Appendix A summarizes necessary facts on the rigged configurations and their 
bijective correspondence with the highest paths. A piecewise linear formula for the 
KKR bijection is available in (A. 5). 

Appendix B contains the proof of Proposition 3.4. For any path not necessarily 
highest, it essentially identifies the two descriptions of the conserved quantity in 
terms of the energy and the configuration. The former is related to the soliton 
content and the latter to the string content. 
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Appendix C is devoted to a proof of Proposition 3.7, which assures the well- 
definedness of the direct/inverse scattering map <I> (3.15). The key is to investigate 
the relation between the two rigged configurations corresponding to the paths q^r 
and r g with q and r both being highest. 

Appendix D is the proof of Theorem 3.12, which asserts the hnearization of the 
time evolution in terms of the angle variables. Our main idea is to realize non 
highest periodic paths as a segment of a large highest path having the structure 
p p (S) p (S) • • • and apply the KKR theory to the latter. 



2. Periodic box-ball system 

2.1. Crystals and combinatorial R. We recapitulate the basic facts in the crys- 
tal basis theory [22, 8, 3, 9, 18]. Let Bi the the crystal of the l-iold symmetric 
tensor representation of Uq{A^j^'^). As a set it is given by Bi = {x — {xi,X2) G 
(Z>o)^ \ xi + X2 = I}- The element {xi,X2) will also be expressed as the length 
I row shape semistandard tableau containing the letter i xt times. For example, 
Bi = {|T|,[2]}, ^2 = {[ir|,[l2],[22]}. (We omit the frames of the tableaux here- 



after.) The action of Kashiwara operators fi,ei : B B U {0} {i = 0,1) reads 
{fix)j = Xj — Sj^i + Sj^i+i and {e.ix)j = Xj + dj^i — where all the indices are in 

Z2, and if the result does not belong to (Z>o)^, it should be understood as 0. The 
classical part of the weight of .t = (.ti,,T2) G Bi iswt(.T) = lAi—X2ai = (xi — .T2)Ai, 
where Ai and ai = 2 Ai are the fundamental weight and the simple root of Ai . 
For any 6 e .B, set 

Si{b) = max{m > | e^^fe ^ 0}, ipi{b) = max{m > | /f 6 7^ 0}. 

By the definition one has ei{x) = Xi+i and (pi{x) = Xi for x = {xi,X2) & Bi. 

For two crystals B and B' , one can define the tensor product B^ B' = {b^b' | 
b € B,b' € B'}. The operators Ci, fi act on B (g) B' by 



(2.1) iiib^b') = 

(2.2) fi{b<E)b') = 



e,h®h' M ^Piib) > ei{b') 

h®e^b' iiipilb) <£i{b'), 

fib®b' li ipi{b) > ei{b') 

b(^fib' i{^i{b)<ei{b'). 



Here 006' and 600 should be understood as 0. The tensor product Bi^®- ■ -^Bi,, is 
obtained by repeating the above rule. The classical part of the weight oi b G B for 
any B = Bi^ ^ ■■■ ^ Bi, is given by wt (6) = ((^i(6) - £i(6))Ai = {eo{b) - (po{b))Ai. 

The crystal Bi admits the affinization Aff(B;). It is the infinite set Aff'(_B;) = 
{C^b \ b e Bi,d e Z} endowed with the crystal structure CiiC^b) = C^^^^'" {e,b), 
fiiC^b) = ('^^^'•"{fib). The parameter ( is called the spectral parameter. The 
isomorphism of the aflane crystal AS{Bi) AS{Bk) ^ Afr(Bfe) AS{Bi) is the 
unique bijection that commutes with Kashiwara operators (up to a constant shift 
of H below) . It is the q = analogue of the quantum R and called the combinatorial 
R. Explicitly it is given by R : C^x (g) Cy ^ c,e+H(x®v)y (g, Qd-H(x®y)x with 

Xi = Xi + Qi{x,y) - Qi-i{x,y), yi = yi + Qi-i{x,y) - Qi{x,y), 
Qi{x,y) = mm{xi+i,yi), 
(2.3) H{x(gy) = -Qo{x,y). 
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Here x ® y :^ y ® x under the isomorphism Bi ® Bk — Bk® Bi. The relation is 
depicted as 



y 



-X or y- 



-y 



y X 
For example Bi ~ Si (g) Bi is listed as follows: 



1 

H=0 



2 

H=0 



l-a+1 a-1 l-a a 



l-a-1 a+1 



H=0 (0<a</) H=-l (0<a<() 

Figure 2.1. Combinatorial R: Bi® Bi B^® Bi 

Let w : (a;i,a;2) i— > {x^,x\) be the involutive Dynkin digram automorphism of 
Bi. The combinatorial R enjoys the symmetry: 

(2.4) {oj®w)R = R{u}®uj) onBi®Bk. 

We write the highest element (/,0) e Bi as Ui. The energy function H (2.3) 
is normalized so as to attain the maximum at H{ui ® Uk) = and ranges over 
— min(Z, fc) < H < on Bi ®i Bk- The combinatorial R satisfies the Yang-Baxter 
relation: 

(2.5) (1 ® R){R ® 1)(1 ®R)^iR(E) 1)(1 (g) R){R (E) 1) 

on AS{Bj) (g) Afr(B;) (g) Afi'(Bfc). Let e(6i ig) • • • ig) 6^) = fofc ig) • • • ig) 6i be the reverse 
ordering of the tensor product for any k. The combinatorial R has the property: 

(2.6) Rq = qR onBi^Bk. 

2.2. Time evolution and conserved quantities. We fix the integer L G Z>i 
corresponding to the system size throughout. Set 

(2.7) V = Bf^, V+ = {pGV\eip = 0}. 

Wc will also write Aff('P) = Aff(Bi)'^^. An clement of V (V+) is called a path 
(highest path). The condition iip = on the path = &i (g • • • (g 6l is equivalent 
to the simple postulate 

(2.8) t{l<i<k\bi = l}>i{l<i<k\bi = 2} for aU 1 < k < L. 

The weight of the path p = h\® ■ ■ ■ is given by wt(p) = wt(6i) H h wt(6i,). 

We write wt(p) > (wt(p) < 0) when it belongs to Z>oAi (Z<oAi). 

The periodic box-ball system is a dynamical system on V equipped with the 
commuting family of time evolutions Ti,T2, . . ., which we shall now introduce. 
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Proposition 2.1. For any path p = hi® ■ ■ ■ ®bL G V and I e Z>i, there exists an 
element vi G Bi such that (e = —di — ■ ■ ■ — dL) 

(2.9) C% {C°h ® • • • C°bL) ^ iC^'b'i C^^bi) ® Cvi 

for some C^^b'-^ ® ■ ■ ■ ® C^^b'j^ ^ Aff(7') under the isomorphism AS{Bi) (g) Aff(7') ^ 
Aff('P)ig)Aff(Bi). Such VI is unique except wt{p) = case, where C'^b'i®- • •®C'^^^z, 
is independent of the possibly non unique choice of vi . 

Proof. Suppose the relation {xi , X2)®p — p' ® (j/i, 2/2) holds under the isomorphism 
Bii^V c^V I® Bi. Setting p' = 6'^ (g) • • • (g) 6'^- , we depict it as Figure 2.2. 




Figure 2.2 



Fixing p and /, we regard it as the functional relation 1/2 = V2{x2)- Then from 
Figure 2.1 we see that 2/2(3^2 + 1) = 2/2 (2^2) or 2/2(2:2 + 1) = 2/2(2^2) + 1- Moreover, 
the latter holds only if all the intermediate vertices in Figure 2.2 are the bottom 
two types in Figure 2.1, namely wt(p) = — wt(p'). 

Let M be the number of 2 € Bi contained in p. We first consider the case 
L — M> M, i.e., wt{p) > 0. Set c = 2/2(0). Suppose that 2/2(0) = c and 2/2(0+ 1) = 
c + 1 hold. From the above observation, the number of 2 contained in p' (for 
{xi,X2) = {I — a, a) in Figure 2.2) must be L — M. Then the weight conservation 
demands that a + M = c + L — M. By the assumption M < L — M, this can not 
happen for < a < c. This implies that 2/2(0) = c for < o < min(c + 1,1), which 
confirms the sought assertion. The case L — M < M can be shown similarly by 
interchanging the role of the letters 1 and 2. Next we consider the case L — M = M, 
i.e., wt{p) = 0. By the same argument as above, we find that 2/2(0) = c for < o < c 
and 272(0) = a can happen for some interval c < a < c'(< I). When c < a < c', 
all the intermediate vertices in Figure 2.2 are the bottom two types in Figure 2.1, 
whose energy H and the vertical Bi part are independent of a. This verifies the 
existence of vi and uniqueness of C'^&'i (g • • • (g C'^^b'j^. □ 

We define the time evolution Ti and the energy Ei of a path p by Ti{p) = 
6'i<g- • -^b'j^ (g V) and Ei{p) = e (€ Z>o) using the notation in (2.9). Proposition 2.1 
assures that they are solely determined from p and I. The definition is summarized 
by the relation 

(2.10) Q°vi®pc^Ti{p)®C^'^P'^vi 

omitting the spectral parameters attached to p and Ti{p). Here vi G Bi can be 
constructed from p hy 

ui (E) p — p* (E) Vl if wt(p) > 0, 
^ ^ uj{ui) (g) p ~ p* (g) Vl ifwt(p)<0 

for some p* gV under the isomorphism Bi®V ^^V ® Bi, where = (^, 0) G Bi as 
defined after (2.4). One may either use the latter relation to define vi when wt(p) = 
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0. Clearly the time evolutions are weight preserving, i.e., wt(T;(p)) = wt{p). They 
are all invertible. By using (2.6), the inverse can be found by 

(2.12) Tf\p) = QTi{e{p)). 

In the rest of the paper wc will use, often without explicitly mentioning, the fact 
that (2.10) determines Ti{p) and Ei{p) unambiguously even though vi is not unique 
in general. If p is a highest path, the affinization of (2.11) is available in Lemma 
B.l. 

Note from Figure 2.1 that the combinatorial Ron B-i_® Bi is the identity map. 
It follows that Ti acts as the cyclic shift or 'exp (\/^ (momentum))': 

(2.13) Ti(6i O 62 • • • O &l) = &i O 6i O • • • O 

Ti will play a special role in our inverse scattering formalism in Section 3. 

Theorem 2.2. The commutativity TiTk{p) = TkTi{p) and the conservation Ei{Tk{p)) 
Ei{p) hold. 

The commutativity with Ti (2.13) is the origin of the adjective "periodic". 

Proof. Take Vk for p and for Tfc(p) as in (2.11). Set R{C°vi®C"vk) = C^Vk®C~^vi 
and regard p as an element of Aff (P). By using the combinatorial R, one can reorder 
C^vi (g) (^Vk (Sip in two ways along the isomorphism AS{Bi) (g) AS{Bk) (S) Aff('P) ~ 
AS{V) Afr(Bfe) O AS{Bi) as follows: 




Tlip) 



Tk(j>) 



TkTiip) 




TiTkip) 



where the equality of the result is due to the Yang-Baxter equation (2.5). Here we 
have identified the outputs with TkTi{p) , C.^'^^'^^'^p^^^^ v^, etc. In particular Proposi- 
tion 2.1 guarantees that ® Ti{p) ~ TkTi{p) i^Vk and zJ; (g)p ~ Ti{p)®vi up to the 
spectral parameter. The sought relations TiTi~{p) = TkTi{p) and Ei{Tk{p)) = Ei{p) 
are obtained by comparing the two sides. □ 

2.3. Extended aflflne Weyl group invariance. Let Sj {i = 0, 1) be the Weyl 
group operator [22] acting on any crystal B as 

j^^e.(.)-..(.)(,) ^.(6)<e,:(6) 

for b G B. Wc extend uj introduced around (2.4) to any B = i?;^ (8) • • • (8) by 
u>{B) = u){Bi^)®- ■ ■®uj{Bi^). Then Vr(A^^'') = (w, sq, Si) acts on V as the extended 
afhne Weyl group of type ^4^^'. 



The action of fi,ei and Si is determined in principle by (2.1) and (2.2). Here 
we explain the signature rule to find the action on any Bi-^ (g) • • • (g) Bi^ which is of 
great practical use. It will be the basic ingredient in proving Proposition 2.3. The 

i-signature of an element 6 € .B; is the symbol 1 1-. The i-signature of 

the tensor product bi (g) ■ ■ ■ (g) bk E Bi^ ® ■ ■ ■ ® Bi^, is the array of the i-signature of 
each bj . Here is an example from Br^ (g) B2 ^ Bi (g) B4: 

11112 O 12 O 2 (g) 1122 

0— sisnaturc ^ ^ "t" I ^ 



1— signature 



-++++ -+ ++ 



where 1122 for example represents 1 1122 | G B4, and not [T|(8)[T|(8)[T|(8)[T| e etc. 



In the i-signature, one eliminates the neighboring pair H — (not — h) successively 

a 13 

to finally reach the pattern 1 h called reduced i-signature. The result is 

independent of the order of the eliminations when it can be done simultaneously in 
more than one places. The reduced i-signaturc tells that £,(6i • • • (S) 6fc) = cx and 
fiipi ® ■ ■ ■ ® bj.) = (3 . In the above example, we get 

11112 ® 12®2® 1122 

0— signature ^ 
1 — signature ~^ ~' 

Thus £0 = 4, tpo = 2, £1 = 1 and = 3. Finally hits the component that is 
responsible for the leftmost + in the reduced i-signature making it — . Similarly, 

fii hits the component corresponding to the rightmost — in the reduced i-signature 
making it +. If there is no such + or — to hit, the result of the action is 0. Si acts 

a 13 Pa 

so as to change the reduced i-signature — • • h • — h into — • • 1 — • +. In the 

above example, we have 



p = 


11112$ 


5 12$ 


52$ 


5 1122 


foip) = 


11112$ 


5 12$ 


52$ 


5 1112 


Mp) = 


11122$ 


5 12$ 


52$ 


5 1122 


eo{p) = 


11122$ 


5 12$ 


52$ 


5 1122 


ei{p) = 


11111$ 


512$ 


52$ 


5 1122 


soip) = 


11222$ 


512$ 


52$ 


5 1122 


si{p) = 


11122$ 


512$ 


52$ 


5 1222. 



For both i — and 1, note that wt(si(p)) = — wt(p) for any p, and Si{p) = p if 
wt(p) = 0. In order that eip = to hold for a path p G P, it is necessary and 

sufficient that the reduced 1-signature to become H h, which is equivalent to the 

condition (2.8). Wc note that elimination of the H — pairs is described by successive 
applications of the rule 

a f3 (J3-oc) + 

where (c)+ = max(c, 0). The time evolutions Ti and the energy Ei enjoy the 
extended affine Weyl group symmetry. 



Proposition 2.3. For any w G p £ V and I G Z>i, the commutativity 

wTi{p) = Ti{w{p)) and the invariance = Ei{p) are valid. 

This property persists in the most general An^ case as announced in [2] . 



Proof. Let the reduced z-signature of p be — • • h • — h and the one for T; (p) be 

■ +. Since T; is weight preserving, one has a — /? = a' — /?'. First we 
show the assertion for w = Si. We assume a < f3, hence Si{p) — fi^°'{p) and 
Si{Ti{p)) = f^~"{Ti{p)). The proof for the case a > P is parallel. Let vi G Bi 
be the element specified by (2.10) from p and set a = Si{vi) and b = ipi{vi). The 
i-signature of (2.10) reads 



(2.14) + . 

where | signifies the position of ® that separates vi with p or T; {p) . 

(i) Case a>h. Comparing the reduced + signature, we see (a <)/3 = 6 + (/?' — 
a)+ < a + (/?' — a)+, compelling /?' > a. Thus we find (3' — a> (3 — a. By taking 
these facts into account, (2.14) is reduced to 

a a — b [3 a fj' — a 

In view of this and 0' — a > — a, application of /f ~" to (2.10) yields vi 

/f-"(p) ~ ffiTiip)) ® C""'^"^";- Therefore vi ® Si(p) ~ Si(Ti(p)) ® C'^'^^^^^i, 
saying Ti{si(p)) = Si{Ti{p)) and Ei{si{p)) = Ei{p). 

(ii) Case a < h. Comparing the reduced + signature in (2.14), we see (6 < 
)h + (3-a = h+{P' - a)+, compclhng /3' - a > 0. Thus we find f3' - a = (3 - a. By 
taking these facts into account, (2.14) is reduced to 



a' [i—a 



In view of this, the application of f^'^^ ^" to (2.10) yields C,^v[ ® /f "(p) ~ 
fi^^iTiip)) (8) C^+^'^P'>v'i with v'l = f^-^vi and 6 = {b - a)Si^o. This is equiva- 
lent to v'l (g) Si{p) ~ Si{Ti{p)) ('^''^P^v'i saying again that T;(s,(p)) = s^{Ti{p)) and 

i;i(s,(p)) = £;,(p). 

Next we show the assertion for w = to. Due to the symmetry (2.4), the relation 
(2.10) implies uj{vi) w{p) ~ u){Ti{p)) w(?;0 at least under Bi igi V V Bi 
forgetting the spectral parameter. Thus we obtain Ti{uj{p)) = Lu{Ti{p)). In (2.10), 
imagine the process of sending vi to the right through p along Figure 2.2 using the 
local rule listed in Figure 2.1. By the definition, —Ei{p) is the sum of H attached 
to the intermediate vertices. Upon application of lo, the four types of vertices in 
Figure 2.1 are interchanged horizontally. But this does not alter the sum of H since 
the number of bottom two types are equal because of wt(p) = wt{Ti{p)). □ 
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(2.15) T^p) = 



Example 2.4. The following is a commutative diagram among the paths mV = 
212221 A 211121 A 222121 A 111212 

^2 J, J, J, J, 

121222 A 121112 A 122212 A 211121 

For all the paths here, we have Ei{p) = 2 for any Z > 1. 

The time evolution T; defined by (2.10) has a simple description for I sufficiently 
large. 

Proposition 2.5. For any path p G V, there exists k G Z>i such that Ti{p) is 
independent of I for I > k. Denoting it by T^{p), one has 

^^{soip)) «/wt(p)>0, 
^ijisiip)) z/wt(p)<0. 

In particular, Too{p) = u){p) is valid ifwt{p) = 0. 

The combination ujSi represents a translation in VF(A^"'^''). On account of (3.14) 
and Theorem 3.12, the minimum of such k is given hy k = js with jg specified in 
(3.8). This is the amplitude of the largest soliton involved in p. The key to the 
connection of the combinatorial R and the signature rule is 

Lemma 2.6. Let {I — a,a) ® h ~ h' ® {I — a', a') under the combinatorial R : 
Bi® Bi~ Bi® Bi. IfO<a<l, then a' and b' are expressed as 

(2.16) a' = {a-so{b))+ + Mb), b' = ^{4'°^'^-"^^ (b)). 

Proof. In Figure 2.1, the top right pattern does not occur. Then (2.16) is directly 
checked case by case. □ 

Proof of Proposition 2.5. In view oiuj{si{p)) = sq{lu{p)) for any p, the two formulas 
in (2.15) are equivalent. Henceforth we assume wt{p) > 0, i.e., eo{p) > ipo{p)- 
Notice that a' in (2.16) is the number of + in the reduced 0-signature of {I— a, a)(Sib G 
Bi® Bi. Consider Figure 2.2 with p = bi® ■ ■ ■ ®bL and (xi, X2) = (I — X2,X2) G Bi 
with any fixed p and X2. It determines b'i,...,b'^ and {y 1,112) successively from the 
left by the local rule specified in Lemma 2.6 by taking I sufficiently large. From 
the above fact and the signature rule, we find that the choice (xi, X2) = (/, 0) leads 
to (2/1,2/2) = {I — fo{p)j'fio{p)) which is called Vi in (2.11). Moreover the latter 
formula in (2.16) implies that the choice (xi, 3:2) = (l — a, a) yield b[® ■ ■■ ®b'j^ = 
u){eQ°^^^ °(p)) provided a < eo{p). By the assumption we are allowed to take 
« = '^o(p)) leading to b'^ ® ■ ■ ■ ® b'j^ = u){so{p)). But this must be Ti{p) due to 
(2.10). □ 

Example 2.7. Take the path p = 121221221111 e Bf^^ satisfying wt(p) > 0. To 
compute T(x)(p) = w{so{p)), we display the 0-signature: 

p=l 2122122111 1, 

- (+ -) (+ (+ -) (+ (+ -) -) -) - 

where the (H — ) pairs to be eliminated successively to get the reduced 0-signature 
are indicated by parentheses. Those paired remain unchanged under so, and we 
thus find 

so{p) = 221 22122111 2. 
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Interchanging 1 and 2 here we obtain Too{p) = 112112112221. In this example the 
composition losq has the effect of doing nothing for the unpaired 1 and interchanging 
the paired 1 and 2. When there remain unpaired + in the reduced 0-signature like 
p = 11122, the effect of wsq is described in the same manner if those + are paired 
with — cychcally. Thus the formula (2.15) reproduces the description of Too in 
terms of the "arc rule" in [15]. 

3. Inverse scattering method 

Here we use the rigged configurations and their bijective correspondence (j) with 
the highest paths 7^+ [5, 6] summarized in Appendix A. 

3.1. Action variable. The set of states V = Bf^ of the periodic box-ball system 
is decomposed into the disjoint union according to the value of the conserved quan- 
tities {El I I G Z>i} which wc called energy. Our aim here is to determine their 
spectrum by making a connection with rigged configurations. It will be attained in 
Proposition 3.4 and (3.5). 

Lemma 3.1. For any path p gV with wt(p) > 0, there is an integer d G Z and a 
highest path p+ G V+ such that p = Tf{p+). 

The pair {d,p+) is not unique in general even if d is restricted to < d < L. 

The proof is elementary, and is illustrated along 

Example 3.2. Take p = 2211221112122111221 with length L 19. Regard the 
letters 1 and 2 in the doubled path p p as the arrows y and \ respectively, 
and construct a length 2L trail by following them. In the example, one has the 
following: 




Find a local minimum that are not higher than any other points on its right. In 
view of wt(p) > 0, such minimums can be found within the first L steps of the trail. 
To obtain p+, read the trail from any one of the minimums to the right for L steps. 
In the example, there are three such minimums and accordingly p — T{^(pi) = 
Tf{p2) = Tl^iPa) with pi = 1122111212211122122, p2 = 1112122111221221122 
and P3 = 1112212211221112122. 

Now we introduce 

(3.1) M = {m = {mj)j>i \ mj G Z>o, ^jm-j < L/2}. 

By the definition ruj = for j » 1, and we identify (mi, . . . , m;, 0, 0, . . .) with 

(mi, . . . ,mi). is a finite set. We call its elements action variables, m, = {m,j) 
is identified with the Young diagram in which nij is the number of length j rows. 
We call the mj x j rectangle consisting of the length j rows a block. 
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Based on the KKR bijection we define the map 
(3.2) 

p I — > m 

as follows. Find {d,p+) e Z x such that p = T^{p+) if wt(p) > and u){p) = 
Tf{p+) if wt(p) < according to Lemma 3.1. Then m is obtained as the configura- 
tion part of the rigged configuration (m, J) = 4'{p+)- As noted in Lemma 3.1, the 
choice of {d,p+) is not unique. Therefore to make sense of the above definition, we 
have to guarantee that if Tf(p+) = Tf (p^) and (f>{p+) = {m, J), (f>{p'^) = {m' , J'), 
then m = m' holds. To sec this, note from the highest condition (2.8) that the sit- 
uation T^~'^ {p+) = p'+ can happen only when p+ = r iSi q and p'_^_ = q®r for some 
shorter highest paths q e Bf'^^'^' and r G B®L-d+d' _ (Without loss of generality, 
Q < d — d' < L may be assumed.) Then the assertion m = m' is included in Lemma 
C.3. 

Example 3.3. Consider p in Example 3.2, which can be expressed by the high- 
est paths Pi,P2 and p^- Computing 4>{pi), one finds that they all lead to m = 
(mi, m2, ms) = (2, 2, 1). The result is depicted as follows: 



2211221112122111221 



Let us relate the action variables m, namely the configurations in the KKR 
theory, to the energy Ei{p) of a path p determined by (2.10). 

Proposition 3.4. For any path p gV, its energy is expressed as 

(3.3) Ei{p) =Y,^Hhk)mk 

k>l 

in terms of the action variable m = {ruj) = ii{p) € M. 

The proof reduces to the highest case p e Vj^, and is given in Appendix B. 
According to Remark B.2, m — {nij) has the meaning of soliton content, i.e., there 
are rUj solitons with length j. See also Example 4.5. The right hand side of (3.3) 
is the number of boxes in the first I columns of the Young diagram corresponding 
to m. Proposition 3.4 determines the range of the energy, i.e., the spectrum of the 
periodic box-ball system, in terms of the action variable m & Ad (3.1). In particular 
it implies the property: 

< E^{p) < E^ip) <---<E,{p) = Es+i{p) = • • • = M, 

where M = J2k>i * greatest integer such that nig > 0, or equiva- 

lently, the length of the longest (top) row of the Young diagram for m. The relation 
(3.3) also determines m from {Ei}. Combining Proposition 3.4 with Theorem 2.2 

and Proposition 2.3, we obtain 

CoroUctry 3.5. The action variable is invariant under the time evolutions and the 
extended affine Weyl group. Namely, ii{Ti{p)) = ii{w{p)) = nip) for any I G Z>i 
andw €W{A^^^). 
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Vim) 



For each m = (rrij) G A^, we introduce the corresponding "level set", namely, 
the set of paths p characterized by = m. Rephrasing the condition in terms of 
the energy by Proposition 3.4, we put 

Vim) = {pGV\Eiip) = ^ min(Z, k)mk for any I}, P+(m) = Vim) n V+. 

k>l 

One has the disjoint union decomposition: 

(3.5) V = U„,eMVim), V+ = U^eMV+im). 

From Corollary 3.5, we see that each Vim,) is invariant under time evolutions 
Ti as well as the extended afhne Wcyl group W''(aJ^-'). (On the other hand highest 
paths can not remain highest under them in genca'al.) In the KKR algorithm, 
M = X^fe>i kmk is the number of 2 G i?i contained in a highest (hence wt(p) > 0) 
path. Thus we find 

wt(P(m)) = {pooAi, -PooAi}, wt(P+(m)) = {pooAi}, 

where = L — 2M > is the limiting (minimum) value of the vacancy number 
(3.7). Thus we have a further decomposition with respect to the weights: 

>(m)Uw(P(m)) ifL>2M, 
Vim) if L = 2M, 

where M = X]fc>i kmk as above and 

Vim) = {pe Vim) I wt(p) >0} = {pe Vim) \ wt(p) =PooAi} 

is a fixed weight subset of Vim). One has Vim) D V+im) and 

(3.6) {p&V\ wt(p) > 0} = UmeMVim), 

\Vim)\ = l^i-)' 
|wt(P(m))| 

for any m e M.. This cardinality will be evaluated explicitly in (4.21) and (4.8). 
The set Vim) is still invariant under any time evolution T;. Now Lemma 3.1 is 

refined into 

Lemma 3.6. For any path p € Vim), there is an integer d G Z and a highest path 
p+ G V+im) such that p = Tfip+). Conversely, for any p+ G V+im) and d G Z, 

one has Tfip+) G Vim). 

3.2. Angle variable. Here we construct the set of angle variables Jim) for each 
prescribed value of the action variable m G M.. We introduce the vacancy numbers 

(3.7) pj = min(j, k)mk. 

k>l 

Given an element m. ^ M, put 

(3.8) H := {j G Z>i I m, > 0} = {ji < j2 < • • • < Js}- 

The vacancy numbers satisfy L = po > Pi > P2 > • • • > ft, = • • • = Poo > 0. We 
introduce 

(3.9) J = Jim) = Jj, xjj,x---x Jj^ , 

(3.10) Jj = {(Ji)iez I G Z, Ji< Jj+i, Jj+„^. = Ji+Pj for all i}. 
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Note that only nij of (Jj), say Ji, J2, ■ ■ ■ , Jruj, are independent, and the Young 
diagram for the partition {Ji+mj — Ji, Ji+mj-i — Ji, ■ ■ ■ , Ji+i — Ji) is contained in 
the nij X pj rectangle for any i. 
For k we introduce a map 

(3.11) (Ji)iez ^ {JDiez, 

Ji = Ji+Sj,k +2min(j,fc). 

We extend ak to the map J ^ J via Ok{J) = crkiJji) x crkiJj2) x • • • x ak{JjJ- 
Obviously the inverse ct^^ exists and akcri = aiCJk holds for any k,l G Z>i. Thus 
the set of elements {crj^f^j^ ' ' " (^j'^ I 'ni,n2, ■ ■ ■ ,ns S Z} forms an infinite Abelian 
group A isomorphic to Z*. We call an element of ^ a slide, having in mind the 
index shift i i + Sj^k in (3.11). The origin of this curious map will bo clarified 
in Lemma 4.2 in connection with the Bethe ansatz. Given two elements J,K & J, 
say that J and K are equivalent and denote by J ~ ii' if J = (y{K) for some a €: A. 
Now we define our main object J = J{m) by 

(3.12) J = J/c^. 

An element of J is called an angle variable. For J G J, we will mostly use the 

same symbol J G ^7 to denote its image in and [J] when emphasis is favorable. 
They are the data of the form ,/ = (Jj^"*^) with i G Z and j G H. For J e J (resp. 
J G »/), it is easily seen that J + d := + (i)igz with d G Z remains in J (resp. 
J). For J G J{m) it is a good exercise to show 

(3.13) al^^---aJ;^{J) = J + L. 

Therefore [J] = [ J + L] represent the same angle variable. JT" is a finite set. Its 

cardinality, namely the 'volume of the iso-level subset in the phase space' will be 
determined in (4.21) and (4.6). We introduce the time evolution of angle variables 
as follows: 



Ti : J{m) J{m) 
^^'^^^ {J^%ezjeH ^ (ji^'^+min(i,Z)),ez,,eH. 



This is a linear flow on the angle variables. It enjoys the commutativity TiTk = TkTi. 
The T; here will be identified with the time evolution (2.10) on the paths V in 
Theorem 3.12. Note that Tf{J) = J + d hence T/^( J) = J G J{m). We also define 
Ti on ^{m) by the same formula as (3.14). 

It is convenient to depict the element {jl''^)i^z,jeH G ^{m) as the Young dia- 
gram m whose rows are assigned with these numbers as follows: 
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This finite data is enough to recover the whole sequence ( J, 



0> 



on account of the 



quasi-periodicity Ji^ruj ~ '^i'" '^Pj (3-10). Then the sUde crj^ (3.11) translates the 

infinite sequence {Ji^''^)i^z by / beside the block-wise shift 2mui{jk,a) on J^. It 
is helpful to depict the action as 



3k 



jUs) 



3k 



2jk 



jOfc-i) 



Here A, 



jUi) 



^Uk - 3)1 



jUa ) 



jUk+i) 



Observe that a\ 

Jk 



only causes the uniform shift 2jkl on 



the upper blocks Jji^.^.^ , ■ ■ ■ , Jj, , whereas it relatively induces the back flow Aq, on 
the lower blocks J^a = >Jji, ■ ■ ■ ,>Jjk-i- ^ rigged configuration {m, J) with J = 
i'^i''^)i<i<mj,jeH G Rig(m) (A. 3) also has the data structure that can be depicted 
as in the above picture. 

3.3. Direct and inverse scattering transforms. The direct and inverse scatter- 
ing transforms are the maps between the path V and the action-angle variables. In 
Section 3.1 we have constructed the map /x (3.2) to find the action variable m M. 
for a given path p G V. It can be found either from the energy Ei{p) by (3.3) or 
from the rigged configuration of a highest path described under (3.2). It remains 
invariant under any time evolution as noted in Corollary 3.5. 

In the remainder of Section 3, we shall only consider angle variables supposing 
that the action variable m has been determined. Wo assume that wt{p) > and for- 
mulate the inverse scattering method for each ■P(to) appearing in the decomposition 
(3.6). The other case wt{p) < is reduced to it by 

Ti{p)=ujTiiLo{p)), 



owing to Proposition 2.3. 

First we formulate the direct scattering transform $ as follows: 



(3.15) 



r{m) 
p 



Z X V+{m) 
{d,P+) 



J{m) 
i{J) + d 



J{m) 
[i{J)+d\ 



Here the pair (d, p+) is the one satisfying p = Tf{p+) whose existence is assured 
in Lemma 3.6. Then the rigging J e Rig(m) is specified by the KKR bijection 
0(p+) = (m, J). The appearance of m here is guaranteed by Lemma 3.6. The map 
L is defined by 



(3.16) 



)l<i<m 



Rig(m) 



J{m) 



l<i<m 
15 



,0 - ((# 



where the infinite sequence {J^ is the one that extends {Ji )i<i<mj quasi- 

periodically as = jf"*^ + Pj for all i € Z. The resulting sequence {Ji''^)iei 

automatically satisfies (3.10). The map t is just the embedding of the quasi-periodic 
extension, hence an injection. 

Note that the {d,p+) satisfying^ = Tf{p+) is not unique for a given p. Therefore 
to assure the well-definedness of one has to show the => part of 

Proposition 3.7. Let p+,p'^ E V+{in) be the highest paths and J, J' G Rig(m) be 
the corresponding rigging, namely, (t>(p+) = {m, J) and 4'(p'+) = ("T'i J')- Then the 
following relation is valid: 

(3.17) Tf{p+) = Tfip'^) ^ i{J) + ~ l{J') +d'€ J{m), 

The proof is available in Appendix C. 

Example 3.8. Consider the length L = 25 path p = 2122112211221111222111122. 
According to Lemma 3.1, it can be expressed as p = Tl'^{p^) = Tl^{p'_^) in terms 
of the highest paths: 

p+ = 1111222111122212211221122 p+ = 1111222122112211221111222. 

One applies the KKR bijection to get </>(?>+) = (m, J) and (f>{p'^) = (m, J'). The 
elements i.{J) + 12 and l{J') + 19 of i7(m) are depicted as 

Tl o 

_zr3 

_ 4 
_ 4 
5 

where 1, 3, 7, 15 on the left are the vacancy numbers (3.7) exhibited for convenience. 
The equivalence l{J) + 12 ~ l{J') + 19 is realized by a3{i{J) + 12) = t(J') + 19. 

Next we show that the direct scattering map $ (3.15) is invcrtiblc, which yield 
the inverse scattering map For this two properties are to be established. 

First, any element in ^{m) must be equivalent to the form l{J) + d for some 
J € Rig(TO) and d E "L, which ensures the existence of an inverse image for the 
middle arrow in (3.15). (</>(?<+) = (m, J).) Second, any two equivalent forms 
l{J) + d ~ I'iJ') + d' G J{m) must be pulled back to the same path in 'P(m) in 
(3.15). The second property is nothing but the part of Proposition 3.7. Thus it 
remains to verify the first property. This is done in 

Lemma 3.9. For any J G J{m), there exist d G Z and J S Rig(m) such that 
J ~ «,( J) + d. 

Proof. Let H = {ji < ■ ■ ■ < jg} be the list of lengths of rows in m as in (3.8). 

We give a concrete algorithm to repair J by a slide so that a"^ ■ ■ ■ a"' J = l{ J) + d 
holds for some J G Rig(m) and d £ Z. Such a slide is not unique in general 
but this relation can always be achieved gradually from the longer rows in m as 

explained below. First concentrate only on the two blocks J^j^_i and J^j_^. In the 
slide crj3_j^ (3.11), the quantity 2min(j, js_i) is a common constant 2js-i on them. 

Therefore by applying cr"^^^ for some n, one can make the resulting ( Jj''"'°"^')i G 

Jj-^_, and (J|-'''^)i G Jj, to satisfy d < J^^'-'^ < ■■ ■ < Jm,;'l < d + Pj^_, and 

d < Ji < ■ ■ ■ < Jmji < d + pj^ for some d. The effect of the common change 
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12 + 



1 
3 

7 

15 



10 



19 



1 
3 

7 

15 



2niin(j, can be absorbed into d. Next one uses aj^_2 similarly to adjust J'j,_2 
to the amended jT^-^^i , J}, up to a redefinition of d without violating the foregoing 
adjustment among the latter two. This process, although the adjustment is not 



unique in general, works through until d < J{ 
for all j G H for some d. 

Example 3.10. 

1 [ 
3- 



U) 



< 



< M < d ■ 



■ Pj is achieved 
□ 



5 
5 
10 
14 



6 
6 

-9 
-13 



~ -9 + 













3 






3 




-5 




-9 



~ -3 



Here the vacancy numbers 1, 3, 9 have been shown only in the leftmost diagram. 

From Lemma 3.9 and Proposition 3.7, it follows that the $ in (3.15) is a well- 
defined and invertible map. Taking the disjoint union over the action variable 

mdz A4, we obtain 

Theorem 3.11. The map $ in (3.15) gives the bijection among the set of paths 
LimeM'Pifn) and the set of action-angle variables Umg^Km, J) \ J G J^{m)}. 

3.4. Solution of initial value problem. Now we present our main theorem in 
this paper. 

Theorem 3.12. The following commutative diagram is valid: 



(3.18) 



r{m) 

Ti 

V{m) 



J{m) 

Ti 

J{m) 



Here Ti on the left and the right are given by (2.10) and (3.14), respectively. 

The proof is included in Appendix D. Since the map $ is invertible. Theorem 3.12 

completes the solution of the initial value problem of the periodic box-ball system 
by the inverse scattering method. The time evolution on the paths P{m) has been 
linearized in terms of the angle variables J{m). The number of computational 
steps required for executing <I>~"'^ oTf o <^ is independent of t. 

It also contains the solution of the initial value problem in the box-ball system on 
the semi-infinite lattice Bi • • • as the case L — > oo . This limit is well-defined 
and drastically simplifies our construction so far. For any path &i &2 • • • obeying 
the boundary condition &j = 1 for « :» 1, one can make 1*^" ®b\®b2® ■ ■ ■ highest 
for some n. Thus the degree d e Z in (3.15) is always frozen to = 0. In view of 
Pj oo as L oo, we do not make a non-trivial identification (3.12) under any 
slide (3.11). Consequently (3.15) just becomes $ : p J, which is nothing but 
(j). Namely, the direct and the inverse scattering transforms are the KKR bijection 
itself (without an upper bound on the rigging), reproducing the results in [23, 24] 
essentially. 

Example 3.13. Let us take the length L = 19 path p in Example 3.2 and derive 
the time evolution 

(3.19) Tl°^°{p) = 1211221112122211221, T^°°°{p) = 2112221211221112112 
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based on the inverse scattering formalism in Theorem 3.12. As noted in Exam- 
ple 3.2, we have the expression p = Tl{pi) in terms of the highest path pi. By 
computing the image of the KKR bijection </) oi p\, one finds 



9 — 



1 



The vacancy numbers displayed on the left of the diagram for convenience will be 
omitted in the sequel. By using the linearized time evolution (3.14) one has 



J 2001 
2001 



2000 
1008 
1004 



~ 2446 + 













1 











7 






3 





1 



The last rigged configuration corresponds to the highest pathp' = 1122112112211121222. 

Therefore the map sends the above scattering data to rf'''*^(p') = Tl'^ip') yield- 
ing the first result in (3.19). Similarly the calculation of T^'^^'^{p) goes as follows: 



^1000 



$(p) 



3001 



^740^667 



2001 

2000 
1008 
1004 



7150 + 



The last rigged configuration corresponds to the highest pathp" — 1211221112112211222. 
Therefore the map sends the above scattering data to Tl^^^{p") = Tf{p") 
yielding the second result in (3.19). 



Example 3.14. We take the length L = 26 path p ■■ 
and compute T^^^{p). p is already highest. 



12112211122211121112211111 



10 
16 







390 



265 
260 
138 
130 



^39^26 



572 + 



Since 572 : 



22 X L, we find Ti^°{p) 



p. 



4. Relation with Bethe ansatz at g = 

Here we show that our inverse scattering formalism originates in the Bethc ansatz 
at q = [7]. The angle variables stem from the logarithmic branch of the string 
center equation and the time evolution is the straight motion of its solution. 
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4.1. Bethe ansatz at q = 0. Let us quickly recall the relevant results from 
the Bethe ansatz at g = 0. For the precise definitions and statements, we refer to 
[7]. The Bethe equation for the spin 1/2 one dimensional XXZ chain on length L 
periodic lattice reads [10] 

/ sin7r(Mi + \ -pj sin7r(ui — Uj + 2y/—llij 

ysin7r(Mj — v^^fi) J sin7r(ui — uj — 

for 1 < i < M, where < M < L/2 is the number of down spins preserved by the 
Hamiltonian. The system is associated with the quantum affine algebra Uq{A^^) 
with q = e~^'^^. Fix m = {nij) G A4. String solutions are the ones in which 
{ui, . . . , um} are arranged as 

(4-2) U U [j {u^J^ + V^{j + l-2k)h + elil\l<k<j}, 

j>l l<a<mj „0)gR 

where M = J2j j'^^j and ^afe stands for a small deviation. Ua^ is the string center of 
the a th string of length j. In this context, the data m € is referred as the string 
content. Let H be as in (3.8). For the generic string solution, the Bethe equation 
is linearized at g = into a logarithmic form called the string center equation: 



(4.3) Aj„,kpuf'> = -{pj + rrij + 1) mod 

keH i3=i 



for j ^ H and 1 < a < rrij. Here pj is the vacancy number (3.7) and 

(4.4) Aja,k0 = Sj^kSaAPj +mj) + 2 mm{j, k) - 6j,k- 

From (4.11), the matrix A = (Aj^.k/:)) is invcrtiblc midcr the condition m G A4. 

There are a number of conditions which the solutions of the string center equation 
(4.3) are to satisfy or to be identified thereunder. First, the Bethe vector depends 
on Ui only via e'^^V^Ui ^ Therefore the string center should be understood as 
Ua ^ G M/Z rather than R. Second, the original Bethe equation (4.1) is symmetric 
with respect to wi, . . . ,um, but their permutation does not lead to a new Bethe 
vector. Consequently, we should regard 

for each j. Last, we prohibit Ua^ = u"^ for 1 < a ^ P < mj for any j. This is 
a remnant of the well-known constraint on the Bethe roots so that the associated 
Bethe vector does not vanish. To summarize, we consider off-diagonal solutions 
{ua^) to the string center equation (4.3) that live in 

(4.5) («P , ,...,uli])e {{R/Zr^ - 5m, ) l&m, for each j, 

where 5n = {(«!,...,««) S (M/Z)" \ Va = for some 1 < a ^ (3 < n}. For 

simplicity wc will often say Bethe roots to mean the off-diagonal solutions to the 
string center equation. Let U{m) be the set of the Bethe roots having the string 
content m. 
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For m € M, we introduce 

(46) o(".)=we.f)n^(''':;-r')- 

(4.7) F = {Fj^k)j^keH, Fj^k = Sj,kPj + 2 miii(j, k)mk. 

In case H = $ (i.e., m = (0,0, . . .)), we put Q.{m) = 1. By expanding the deter- 
minant, it is easy to sec 0(to) G Z. Moreover, (4.13) below tells that f2(m) € Z>i 
and (4.6) is also expressed as 

(4.8) 0(m) = — n 



m 



where pj, = L — 2M with M = "^jj^nj. In case p^^, = 0, the combination 
J^fPjs+"-i'js-^\ is to be understood as Note that Cl(m) = detF given in 

(4.13) if m^- = 1 for aU j e iJ. 

Theorem 4.1 ([7] Theorems 3.5, 4.9). 

(4.9) n{m) = \U{m)\ mGM, 

(4.10) J2 "("^) = (m) < M < L/2, 

where the sum extends over mi, m2, . . . € Z>o smc/i t/iaf jruj = M. 

The derivation of (4.9) is due to the Mobius inversion trick. Actually (4.10) is 

known to hold for any > if the symbol (") is interpreted as the generalized 
binomial coefficient a{a—l) ■ ■ ■ (a— /?+l)//?!. The combined identity J^m I ^("^)l ~ 
(jj^) is called the combinatorial completeness of the string hypothesis at q = 0. 
We include the formulas needed here and in Sections 4.2 and 4.3. Fix I e Z>i. 



^'°^'"+^ '"^ detF (0<n<s-l). 

Here the matrix A[kP] is obtained from A = {Aja^k0) (4.4) by replacing the kp th 
column by h. In (4.14), we are using the notation io = and ^„ = min(Z,_7„) for 
1 < n < s. Pif, = Pq = L as noted after (3.8). The matrix F[k] is obtained from 
F = {Fj^k) (4.7) by replacing the column foi k £ H by the column vector h'. (We 
understand F[jo] = 0.) Here the vectors h and h' are given by 

(4.15) h - (min(j, l)),eH, i<a<m, e Z™" +•••+'"- , h' = (min(j, 0),eff E Z^ 

which are dependent on I. Note that (4.14) is vanishing unless < n < t, where t 
is the maximum integer such that it+i > it- By the definition < t < s—1. The 

20 



(4.11) 


det^ 


(4.12) 


det A[kP] 


(4.13) 


detF 


(4.14) 




detF[j„+i] 


-deti^[i„] 



relations (4.11) and (4.13) are derived by noting 

(4.16) Yl ^^a,k0 = E ^^'i' = ^■ 

keH,l<l3<mk keH 

See also eq. (3.8) and (3.10) in [7]. Combining (4.11) and (4.12), we find 

detA[ja] detF[j] 



(4.17) 



{j € H, 1 < a < nij). 



det A det F 

The property detA ^ is obvious from (4.11) and (4.13) under the condition 
m e M. 

4.2. Solitons as strings. Let us uncover the origin of our inverse scattering for- 
malism in the Bethe ansatz. Note that the action variable m G A4 has emerged 
in two independent contexts in J{m) and U{m). In the former it represents the 
soliton content of a path whereas in the latter it is the string content of the Bethe 
equation. The key to their link is the map: 



(4.18) 



J{m) 
J 



U{m) 



where u = (ua^)jeH,i<a<m3 is determined from J = {J!>^')i^z,jeH € J{m) as the 
solution of the linear equation: 



(fe) 



liPj + ruj + 1) + +a-l {jeH, l<a<mj). 



We write it as Au = c+ J + p, where c = (ca'') corresponding to {pj +mj + l)/2 
is a constant vector having the components that are independent of the index a. 

p = {pa^ = a — 1) is also a constant vector^. The equation (4.19) is the string 

center equation (4.3) that corresponds to a prescribed logarithmic branch of the 
Bethe equation. To make sense of (4.18), we are to check the off-diagonal condition 
(4.5) on u. To do this it is useful to grasp the structure of the matrix A. For 
example in case H = {1, 2, 3}, it looks as {Pj = pj + mj for short) 



1 

2 



Vi + l 
2 



2 

^2-1-3 

3 
4 



v 



P2 + 3 4 
4 P3 + 5 



r:i + 5 J 



which consists of 9 sub-matrices of size nii x rrij (1 < i,j < 3). Thus (4.19) leads 



to {pj + mj){u 



J. 



J. 



U) 



a — p. From this and the condition 



(3.10), we confirm that ^ produces a specific array of real numbers {ua ) that obey 



The vector c plays a role only in Proposition 4.11, and a — 1 in (4.19) can be replaced with 
a + a for any integer a. 
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< < a+1 for some a G K, which may indeed be viewed 
as a representative element oi'U{m). 

Varying J G Jim) to J' changes the image u = ^'(J) to u' = 'i'(J') as an 
array of real numbers. But u = and u' = {u'a^) can be regarded as the 

same element in U{m) under the identification by &mj in (4.5) ^. Prom the above 
property of such an event takes place if and only if u' = ■ ■ ■ a^^ {u) for some 
ni,. . . ,ns e where ak is defined by 



)\ ) y—i J • ■ • 7 '^"ij ; + 1) if j — k, 



Since the matrix A is invertible, the effect of dk is translated to that on J e Jijn). 

Lemma 4.2. If u = ^(J), then ai.{u) = ^{ak{J)), where is the slide defined 
in (3.11). 



Proof. Set u' = ak{u) = 4* (J'). By the definition we have Au = c+ .J + p and 
Au' = c+ J' + p. In view of the structure of A, we find for j ^ k that j'a'^^ = Ja ^ + 
^ja,kmk = Ja^ + 2 niin(j, fc). On the other hand, the k th block of c+ J' + p = Au' 
is evaluated as (c = c^^ ) 



f(k) 



c+J[ 

0+4"^ + 1 



Xc+frnJ+mk-lJ 
/Pk+mk + ^k-l 



2k-l 



2k-l 



Pk+rrik + ik-lJ 



,(fe) 



■ E 



Ak2,jf3 



\-^kmk,iP/ 



,0) 



c+J^^^ + l + 2k-l \ 



c+jS+mfe-l + 2fc- 1 
Vc + jf ) + Pfe + TOfe + 2fc - 1/ 



where the last equality is due to the k th block of the relation Au = c+J+p. Noting 

'2k 



the quasi-pcriodicity j}*^-* +pk = J^l+i^ this result is expressed as J\ 
for 1 < i < nik- Unifying the formulas we find that J' = (4''^) is obtained from 



(k) 



f{k) _ jik) 



'U) 



2 min(j, A;) in agreement with (3.11). 



□ 



Wc have established that ^(J) and 'I'(J') represent the same element in lA{m) 
if and only if J ~ J' G J{m). From the definition (3.12), we obtain 



■^A proper treatment of this is to introduce the set of array of real numbers U{m) that projects 
onto U{m) under the identification scheme in an analogous way from J'{m) to J{m). However 
we skip it here supposing no confusion might arise. 
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Theorem 4.3. The map \E' (4-18) induces the bijection between the set of an- 
gle variables J{m) and the set of off-diagonal solutions U{m) to the string center 
equation (4-3). 

The induced bijection will also be denoted by 'I'. Now we are able to introduce 
the time evolution Ti of the Bethe root u = *(J) G U{m) by T;(u) = 4'(T;(J)) 
using (3.14). 

Corollary 4.4. Theorem 3.12 is extended to the following commutative diagram: 

V{m) — ^ J{m) : U{m) 

(4.20) T,| T,| T,| 

V{m) 1 J{m) U{m) 

Their cardinality is given by 

(4.21) |P(m)| = |J(m)| = \U{rn)\ = 0(to). 

The result |7'(?7?,)| = ^{m) endows the character formula in [7] with the quasi- 
particle interpretation. The quantity Q,{m) originally associated with the string 
content m turns out to be the number of states in the periodic box-ball system hav- 
ing the prescribed soliton content m. Then the identity (4.10) implies that weight 
spaces of the s?2 module (C^)®^ are decomposed into subspaces spanned by iso- 
energy states in the periodic box-ball system. These observations and Proposition 
3.4 arc quantitative supports to identify strings with solitons. The identification 
also agrees with the physical picture [12] that the strings represent bound states 
of magnons over the ferromagnetic vacuums 111 . . . Ill or 222 . . . 222. The fact 
[^(m)! = Q.{m) has also been shown in [15] by a different approach. 

Example 4.5. Consider the length L = 8 paths p with wt{p) = 0. There are 
(4) ~ '^^ s^'^'^ paths, which are grouped into T'{m) according to the 5 elements 
m = (mi,m2,...) € M (3.1). Their cardinality |7'(to)| = i}{m) are given as 
follows: 



m= (mi, m2, ma, m4) 


(4,0,0,0) 


(2,1,0,0) 


(0,2,0,0) 


(1,0,1,0) 


(0,0,0,1) 




2 


24 


4 


32 


8 



The paths in V{m) are listed as follows: 



(mi,m2,m3,m4) 


V{m) 


(4,0,0,0) 


12121212, 21212121 


(2,1,0,0) 


Tf (12121122), Tf (12112122), Tf (11212122) 


(0,2,0,0) 


11221122, 21122112, 22112211, 12211221 


(1,0,1,0) 


Ti"(12111222), Ti"(11211222), Ti"(11121222), Ti"(11122122) 


(0,0,0,1) 


Ti"(11112222) 



where Ti{p) stands for the 8 paths {n E Zg) obtained by cyclic shifts of p. The 
path 11122122 for m = (1,0,1,0) is the intermediate stage of the collision of the 
length 3 and 1 solitons. 
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4.3. Periodicity. The time evolution T; is invertible as seen in (2.12) and the set 
of states V is finite. Therefore every path p gV possesses the property T/^ip) = p 
for some integer > 1. We say any such integer a period of p. The minimum 
period is called the fundamental period of p and denoted by Af* = Af^ (p) . Every 
period is a multiple of the fundamental period Af* . Here we establish a formula for 
the fundamental period under any Ti taking advantage of the linearization scheme 
(4.20). We assume p G Umg^7'(m) with no loss of generality thanks to (3.6) and 
Proposition 2.3. 

For p e ■P(m), let J = $(p) e J{m). The fundamental period Af* of p under Ti 
is a function of /, the action variable m and the angle variable J. To be expository, 
we approach the fundamental period in two steps. First wc show in Theorem 4.6 
that there exists a value Af — Afi (m) that is independent of J nonetheless making 
Ti^{p) — p hold for all p e V{m). We call it the generic period inherent to the 
whole V{m). In general Af is yet a multiple of the fundamental period TV*. Second 
we analyze the accidental symmetry gained by special J that makes N* a divisor 
oiAf. It improves the idea of Theorem 4.6 and leads to the final result in Theorem 
4.9. We begin by presenting the formula for the generic period Af. 

where the matrices A, A[ja] and F,F[j] are defined in (4.4), (4.7) and after (4.14). 
In particular, A\ja\ and F\j] are dependent on I. The union (J^.^ extends over 
those j G H, 1 < a < rrij such that det A[ja] ^ 0. The union [j'j does over those 
j E H such that det F[j] ^ 0. The least common multiple LCM(ri, . . . . r„) of 
nonzero rational numbers n , . . . , is the minimum positive rational number in 
Zn n • • • n Zr„. Thus LCM(1, n, . . . , r„) denotes the minimum positive integer in 
Zn n • • • n Zr„. The equality of the two expressions is due to (4.17). 

Theorem 4.6. For any path p G V{m), Tf^{p) = p is valid. 

To be instructive we present two proofs utilizing U (m) and J{m) in the lineariza- 
tion scheme (4.20) although they are essentially the same. In both proofs we first 
consider a sufficient condition for Tf^{p) = p hy regarding Af as an unknown, and 
confirm afterwards that the choice (4.22) fulfills it. The identities and definitions 
(4.11) (4.15) work effectively. 

Proof based onhi{m). Suppose p J ^ u and Tf^{p) ^^> J' u' . We have 
Au = c + J + p and Av! = c + J' + p. From (3.14) and (4.15), this leads to 
A(u' — u) = J' — J = Afh. A sufficient condition for Tf^lp) = p is that all the 

components of v := il' — u arc integers. Since v = AfA~^h, its ja th component is 
v'^i^ = Af dot A[j a]/ det A. Therefore the Af in (4.22) satisfies the condition. □ 

Proof based on J{m). Take J G J{m) such that [J] = $(p) G J{m). We are to 
consider a sufficient condition for [T^{J)] = [J]. This indeed happens if there exist 
integers ui, . . . ,ns such that cr"^^ • • • fT"/(J) = T{^{J) G Jim). We look for them 
by further presuming that ra = Ua/mj^ is an integer. Then from (3.10) and (3.11), 
the above relation on the block Jj^ is expressed as 

s 

(4.23) pji^rk + 2^inin{jk,ja)mj^ra = Af inin{jk,l). 

a=l 
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In terms of the column vectors r = *(ri, . . . ,rs), h' (4.15) and the matrix F (4.7), 
this is written as Ff = J\fh'. Its solution is given as ra = Af det F[ja]/ det F by using 
the matrix F[k] appearing in (4.12). Thus the J\f in (4.22) assures the existence of 
Vra e Z. □ 



The expressions (4.22) can be simplified. We employ the notation used in (4.14), 
e.g., H = {ji, . . . ,js} as in (3.8), «„ = min(j„, ^) (1 < n < s), zq = and <t < s — 
1 is the maximum integer such that it+i > if Since F[jt_^_i] = F[jf^2\ = ■ ■ ■ = F\js] 
from (4.14), we have 



(4.24) 

/ dctF dctF dct^^ 
Af = LCM 1, 



= LCM I 1, 
= LCM 



det F[ji] ' det F[j2] ''"'det F[jt+i] 
det F det F det F 



det F[ji] ' det F[j2] - det F[ji] '"''det - det F[jt] 



where the second equality is an elementary property of LCM. The last line is 
obtained by (4.14). The imion U'^ extends over those n such that is 
finite. This caution is needed only if wt(p) = and I > js, where one encounters 
Pis — Pjs = in the denominator. In this case t = s — 1 and only n = s — 1 is 
allowed in the union so as to cancel the zero by Pi^^.^ = in the numerator. This 
leads to J\f = LCM(1, Noting that = — Pj, = "iijs — js-i)'Tij^, 

we get Af = 2mj,. Thus Tf^ip) = p is certainly valid, for Proposition 2.5 and the 
remark following it imply that T^{p) = p- Another simplification of (4.22) occurs 
at ? = 1, where i„ l(n > 0) hence t = 0. In this case we have Af = L, which is 
again consistent with T/" (p) = p for any path p. 

Being able to simplify (4.22) into (4.24) is gratifying. However as we will see 
in (4.31) in Section 4.5, it is instead the expression (4.22) that makes the generic 
period Af conceptual and elucidates the essence of the game. We note that the 
formula (4.22) is the simplest case of the most general one for An^ conjectured in 
eq.(8) in [2]. The expression (4.24) with / = oo was first obtained in [15]. 



Example 4.7. Take the path p = 21121111221122111111222 of length L = 23. 

It belongs to 'P{m) with m = fi(p) = (mi, m2, ms, 7714) = (1,2,0,1) hence H = 
{1,2,4} and s = 3 in (3.8). The vacancy numbers (3.7) are {pi,P2,P3,P4) = 
(15, 9, 7, 5). For 1 < Z < 4, the data ii, 12, iz and t are listed as follows: 



I 


n 


«2 


is 


t 


1 


1 


1 


1 





2 


1 


2 


2 


1 


3 


1 


2 


3 


2 


4 


1 


2 


4 


2 



25 



By convention we also have io = and po = L = 23. According to the last 
expression in (4.24), Af is calculated as 

1 = 1: LCM(l,po) = LCM(1,23) = 23, 
PiPo \ 115 



2 : LCm(i, 

3 : LCm(i, 

4 : LCm(i, 



=LCm(i,— ,15) =345 
345 135 



P2 

PiPo P2P1 



P3 P3 

PlPO P2P1 P2 

2 



Pa) =LCM(i,— ,— ,9j =3105, 
'2; 



LCM 1,69,27, 



621. 



P4 P4 

Actually these generic period Af coincide with the fundamental period J\f* of p 
under the respective time evolutions Ti. The reason will be explained in (4.27). 

Example 4.8. Take the path p = 2122112211221111222111122 of length L = 25. 
This is the same path as that considered in Example 3.8. It belongs to V{m) 
with m = n{p) = (mi,m2,m3,m4) = (1,2,1,1). The vacancy numbers read 
{pi,P2,P3,P4) = (15, 7, 3, 1). This time we calculate J\f along (4.22). The matrix F 
(4.7) reads 



F 



(p\+ 2toi 
2mi 
2mi 
\ 2toi 



2m2 
P2 + 4m2 

4m2 
4m2 



2m3 

P3 + 6TO3 



2m4 

4:1714 

6 7714 
P4 + 81114 J 



/17 
2 
2 

V2 



4 
15 



2\ 
4 

6 

9/ 



If i = 3 for instance, the matrix F[j] is obtained from F by replacing its j th column 
with /!' = *(1, 2, 3, 3) (4.15). Thus F[l], F[2], F[3] and ^^[4] look as 

/I 



4 2 2 

2 15 4 4 

3 8 9 6 
\3 8 6 9/ 

which lead to 

deti^ 



^17 
2 
2 

V2 



2\ 

4 

6 

9/ 



det det F 



(n 
2 
2 

V2 



dcti^ 



4 

15 



2\ 

4 

6 

9/ 



^17 

2 
2 

V2 



4 

15 



1\ 

2 
3 

3/ 



(1,125, 



875 875 875 \ 
"32"' 157' 157/ ■ 



det F[l\ ' det F[2] ' det F[3] ' det F[4] - 

Taking the LCM of them, we get M = 875. We list the result of such calculations 
for 1 < Z < 4 in the following table. 



I 


LCM 


Af 


1 


1, 


25, 


25, 


25, 


25 


25 


2 


1, 


375 
7 ' 


375 

32 ' 


375 

32 ' 


375 

32 


375 


3 


1, 


125, 


875 
32 ' 


875 
157' 


875 
157 


875 


4 


1, 


375, 


2625 
32 ' 


2625 
157 ' 


875 
344 


2625 



For each I, Tf^{p) = p holds. Actually the Af listed here coincides with the 
fundamental period of p. See (4.27) for the reason. 

As in Examples 4.7 and 4.8, the generic period J\f = Ni{m) (4.22) turns out to 
be the fundamental period Af^ip) for a majority of paths p G V{m). The event 
Afi{m)/Af*{p) G Z>2 indicates some extra symmetry in the path p. The structure 
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of our angle variable provides a lucid picture on the nature of the extra symmetry, 
and improves (4.22) into a formula for the fmidamcntal period. 

Consider the block J — {Ji)i^% G Jj{m) (3.10) for j £ H = {ji, . . . ,js}. Seek 
the maximum positive integer g — gj such that 

(4.25) Pl/JlleZ>i, J.,m, -Ji = ^ for any i € Z, 

9 9 - 9 

which is consistent with the condition Ji+ruj = Ji+Pj in (3.10). Wc say that such J 
has the order g symmetry, which means the finer quasi-periodicity than the original 
one by factor g. The case g = 1 corresponds to the previous treatment. 

Theorem 4.9. For any path p € P(m), let {j'^^^)j^H = ^{p) he the angle variable. 
Suppose J^^^ has the order g^ symmetry. Then the fundamental period ofp is given 
by 



(426) ^.=LCM(i.y;^). 



deti^ \ 

n, 

where U'- is the union over those j such that detF[j] ^ 0. 



Obviously Af > Af* and the equality holds if (but not only if) gj = 1 for all 

jeH. 

Proof. A slight modification of the argument around (4.23) suffices. From the 
definition (4.25) of the symmetry of J, to validate [p] = p or equivalently 
' ' ' '^ll ^"^^ ~ ^® necessary for = ria/ {'iTija/9ja) to be an integer and 

to satisfy 

s 

^r/c + 2^min(jfc, ja)— ra = A^min(jfe,/). 

9jk „=l 9ja 

This is precisely (4.23) with replaced by ra/gj^ and J\f by N. Thus iV must be 
chosen so that ra = Ngj^ det i^[ja]/ deti^ be an integer for all 1 < a < s. The M* 
in (4.26) is the minimum of such N . □ 

Beside the generic case V^j = 1, the next simplest situation is ^gj = g{> 2). 
It only happens when L/^ e Z, and corresponds to the path of the form p = q®^ 
for some q G b®^/^ without any such symmetry. (This statement is justified from 
Lemma C.l.) In this case, M* (4.26) coincides with the generic period N (4.22) for 
q, which is consistent with Ti{p) = Ti{q)®^ . 

Example 4.10. We consider the path p treated in Example 3.14. We know that 
T^^^{p) = p, and actually 130 is the fundamental period under T^. Looking at the 
angle variable J = {J^^\ J^^\ J^^^) = $(p) given there, we find that J^^^ and J^^^ 
possess the order 2 symmetry. Thus (51, 52, ffs) — (2, 2, 1) in the notation in (4.26). 
We have H = {1,2,3}, (mi,m2,m3) = (2,2,1) and {pi,P2,P3) = (16,10,8). Thus 
the matrix F (4.7) reads 

^Pi+2mi 2to2 2m3 
F = I 2mi p2 + 4m2 4m3 

2mi 4m2 ps + 67713 j 
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We consider the time evolution T3, for which the vector h' (4.15) reads h' = *(1, 2, 3). 
Replacing the columns of F with this, we get 

/I 4 2\ /20 1 2\ /20 4 l\ 

F[l] = [2 18 4 , F[2] = 4 2 4 , F[3] = 4 18 2 . 
\3 8 14/ \4 3 14/ \4 8 3/ 

Substituting the determinants (det F, det i^[l], det F[2], det F[3]) = (4160, 80, 288, 704) 
into the formulas for the generic period J\f (4.22) and the fundamental period J\f* 

(4.26) , we find {I = 3) 

x^,.^ 4160 4160 4160\ ,^,,/, _ 130 65\ 
= LCM 1, — — , — — , — — = LCM 1, 52, — , — = 260, 

V ' 80 288 704 / V ' ' 9 11/ 

/ 4160 4160 4160\ / 65 65\ 

M* =LCM(1,- — — ,- — — ,— -) =LCM(1,26,— , — ) = 130. 

V ' 2 X 80 2 X 288 704 y V' ' 9 11/ 

Prom (4.22), (4.25) and (4.26), we find 

(4.27) M = M* if GCD(pj,mj) = 1 for all j e iJ, 

where we employ the convention GCD(0,rnjJ = 1 for the greatest common divisor 
when pj^ = 0. Note that (4.27) is a sufficient but not a necessary condition for 
N = Af*. It explains the reason for Af = Af* in Examples 4.7 and 4.8. For I = 00, 
the fundamental period has also been studied in [15]. 

4.4. Bethe eigenvalue. The time evolution T/ in the periodic box-ball system is 
the q = limit of the row transfer matrix T;(^). Its eigenvalues are given by the 
analytic Bethe ansatz [25, 26]. Let Q{0) = Hftli sinh7r(^— \/^Ufc) be Baxter's Q- 

function, where {uk} satisfy the Bethe equation (4.1). We set q = e~^'^^ as before, 
^ = e^'^^ and assume M < L/2. For the string solution (4.2), the q ^ limit of the 
Bethe eigenvalue ([26] eq.(1.8)) coincides with that of the top term ([26] eq.(2.12)): 

(4-28) i^ S + S ^^"''''^'' A,=exp(27rV^^min(j,0(^^L^V^)) 

under an appropriate normalization. Here Ei = J2k ^^^{h k)mk is the right hand 
side of (3.3), and the sum Y^j^ extends over j G H and 1 < a< rrij. 

Proposition 4.11. For any solution u= {ua^) to the string center equation (4-3), 
the equality = 1 is valid, where Af is the generic period (4-22). 

Proof. Put V = (va^) = AfA~^h, where h is specified in (4.15). Its component is 

expressed as Va'' = A/'dot ^[ja]/ dot A = A/'dct F[j]/ dot _F, where we have used 
(4.17). This is independent of the index a. Moreover from (4.22), we may set 
Va^ = Xj for some Xj G Z. Further define the vectors c = {c^^^)ja with c^^^ = 
{pj + ruj + l)/2, and e = (|)ja- Then the string center equation (4.3) is written 
as Au = / + c for some integer vector /. The A; (4.28) is expressed as A; = 
exp(27r-\/— 1 */i(u + e)) = cxp(27rV— 1 */i^~^(-f + c + Ae)). Since A is a symmetric 
matrix, we have Af{*hA^^)ja — Va^ = Xj. Thus A^ = exp(27rV^^^-^ a;j(/ + 
c+Ae)ja) = exp{2TT\/^YjaXj{c + Ae)ja)- Now (log A^)/(27r\/^) is calculated 

28 



ja fc/3 

^ X] (P^ + rrij + 1 + pj + rrij + 2 ^ niin(j, k)mk - rrij) 

ja k 

a;j?nj' (pj + ^'^^"^ ^ + niin(j, k)mk) = mod Z. 

j k 

□ 

Proposition 4.11 does not serve as a proof of Theorem 4.6. Nevertheless it has 
opened a route to create a conjectural formula for the generic period in a large 
class of generalized periodic box-ball systems [1, 2]. For I = oo, Proposition 4.11 
has been shown also in [27] independently. 

4.5. Discussion. We write the string center equation (4.19) in the matrix form: 

(4.29) Au = c + IeZ^, 

where 7 = mj^ + • • • + ruj^ for H = {ji, . . . ,js} as in (3.8). The shifted array 
I = J + p will also be called the angle variable here. See around (4.19) for p. Prom 
(4.5) and (4.19), the / belongs to the set 

(4.30) {Im^^X---Xl,nJ/r, 

where In = (Z" — A„)/6„ is the n dimensional lattice without the diagonal points 
An = {(-21, ... , Zn) & "L" \ Za = Zfj for some 1 < a ^ (3 < n) identified under the 
permutations S„. F = ®kzH,\<i3<mk '^-^kfi is the 7 dimensional lattice generated 
by the column vectors Akfj of the matrix A = (Aja.kp) (4-4) that characterizes the 
string center equation (4.29)*. The division by F originates in the identification of 
u under u^^^ — > u^^^ + 1 in (4.29). The time evolution T; (3.14) on the angle variable 
/ is expressed as the linear flow Ti{I) = I + h in terms of the Z-dependent vector 
h G IP' defined in (4.15). According to Theorem 4.3, it induces the time evolution 
of the Bethe roots Ti{u) = u + A~^h. They are summarized in the following table: 





Bethe roots u 


Angle variables I 


mod 




r = ..iz" 


Ti 


A-^h 


h 



Here mod is the lattice under which the respective variables are to be identified save 

the permutations &mj^ x ■ ■ ■ x &mj^ ■ Now it is transparent under what condition 
the time evolution Tf^ becomes trivial. According to the above table, it is presented 
in the three equivalent forms: 

(4.31) AfA-^heZ'^, AfheT, Afp-'^h'eZ', 



using (4.4) as 

^^Xj(c + Aejja = 

ja 



■^The column vectors Aj-p are independent because of detj4 > under the condition m & A4 
as noted in the end of Section 4.1. 
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where the last one is a contracted version of the first due to (4.17). See (4.15) for 
h' . Similarly, the formula (4.26) for the fundamental period is rephrased as 

(4.32) M*F-^h' eG-^l^ G = diag(p,-,,...,5,J. 

In this way we arrive at the intrinsic meaning of the generic period (4.22) (rather 

than in the form (4.24)) and the fundamental period (4.26). Namely, they arc the 
smallest positive integers N and N* that make (4.31) and (4.32) valid. Depending 
on the choice of Hn (4.15), the vectors h and h' (like Hamiltonian) encode various 
direction and speed of the straight motions in the set (4.30) corresponding to the 
time evolution Tj. The simplest among them is/i = *(l,l,...,l)€Z''' for Z = 1, for 
which L/i e r holds because of (4.16). This fact corresponds to the simple property 
'^\{p) = P foi' P- See (2.13). More generally, the generic and fundamental 
period imdcr the combined time evolution T = Yli Tf' {Pi G Z) can be obtained 
by replacing h = hi and h! = h[ (4.15) with J2i ^^'^ ^if^i^^i (4.22) and 
(4.26), or equivalently, in (4.31) and (4.32). The set (4.30) is an analogue of the 
Jacobi variety on which the nonlinear dynamics on V looks as a straight motion as 
in the classical theory of quasi-periodic solutions to soliton equations [20, 21]. The 
coefficient A in the string center equation plays the role of the period matrix. Its 
size 7 is the total number of solitons, which is equal to the first energy Ei (3.3). 

Under any set of selected time evolutions Ti^, . . . ,Ti^, one can describe the de- 
composition of V{m) into the disjoint union of orbits: 

V{m) = y {orbit under T,, , . . . , T; J. 

From (4.30), each orbit hero is in one to one correspondence with an element of 

(4.33) (J„^.^ X • • • X )/(r + Zhi, +■■■ + Zhi, ). 

Equivalently, in terms of the original angle variable, the orbits are labeled by 
J{m)/{Tf^ • • • T;^), where the division means the identification of the elements con- 
nected by the time evolution (3.14). 

Under the single time evolution T;, the number of orbits contained in ■p(m) 
is f2(m)/A/' if the action variable m G M satisfies the generic condition (4.27). 
Sec (4.21) and (4.6) for Q{m). Another simple situation is to include the whole 
family Ti,T2, ■ ■ ■■ which maximizes the orbits and minimizes their number. From 
(3.9) (3.12) and (3.14) it is easy to see that J'(m)/(Tfr| • • • ) is factorized as 

(4.34) 4x4x---x4- 

Here Jj = Jjj'^ is obtained from Jj = {(Ji)iez | G Z, Jj < Jj+i, Jj+to^ = 
Ji -Vpj for all i) in (3.10) by the identification (Jj) ^ ( J,-) defined by J[ = Jj+i for 
all i or Jl = Ji + 1 for all i. The difference di = J^+i — Ji + 1 satisfies di £ Z>i and 
dj+i -|- + • • • + di+mj = Pj + fnj, especially di = di+m- for all i £l^. Thus the 
set Jj is in one to one correspondence with the arrangements of pj + rrij letters 1 
and mj letters 2 as 

di d2 ds '^"^3 

that are inequivalent under the periodic boundary condition. No two letters 2 are 
allowed to be adjacent because of rf, > 1. Such arrangements are the states of 

30 



the periodic box-ball system containing amplitude 1 solitons only. In this way we 
obtain 

(4.35) \J{m)/{TfTf ■■■)]= C{pj, , mj,)C{pj, , m,-J • • • C{pj^ , mjj, 

where C(p, m) {p > 0,m > 1) is the number of orbits in the size p + 2m periodic 
box-ball system containing m solitons of amplitude 1 only. Namely, the count of 
orbits under the entire family Ti, ■ • • splits into individual blocks wherein all the 
solitons behave effectively as amplitude 1. The quantity C{p,m) is characterized 
as the number of orbits of the monomials x'^^ ■ ■ ■ a;^™ under the cyclic shift Xi — > 
Xi^i {xi+rn = Xi) in the complete symmetric function hp(xi, . . . ,Xm)- 

Finally we comment on the motion of the Bethe roots w i— > T; (u) = u + A~ ^ h. 
The fusion transfer matrix corresponding to T; does not change the Bethc vectors 
up to an overall scalar nor the associated Bethe roots. This does not contradict the 
motion of u considered here since each path that we associate to u is a monomial 
in (C^)®^, which is not a Bethe vector at g = in general. 

5. Summary 

In this paper we solved the initial value problem in the periodic box-ball system 
by a unification of the combinatorial Bethe ansatze at q = I and q = 0. Section 
2 gives the formulation of the periodic box-ball system in terms of crystal basis 
theory. The commutativity, energy conservation (Theorem 2.2) and the invariance 
under the extended affine Weyl group (Propositions 2.3 and 2.5) in the periodic 
setting are firstly shown explicitly in this paper. In Section 3, we introduced the 
action and angle variables in (3.2) and (3.12) and the linear time evolution (3.14) on 
the latter. The direct/inverse scattering map is defined in (3.15), which linearizes 
the dynamics as in Theorem 3.12. In Section 4, our inverse scattering formalism are 
linked with the Bethe ansatz by the key relation (4.18). The action-angle variables 
are in one to one correspondence with the off-diagonal solutions (4.5) to the string 
center equation (4.3) as summarized in Corollary 4.4. It has led to the explicit 
formula (4.21), (4.6) counting the states characterized either by soliton content 
(energy) or string content (configuration), which we identified in Proposition 3.4 
and Section 4.2. As further applications, the generic (4.22) and fundamental (4.26) 
period and the number of disjoint orbits (4.35) under the commuting family of time 
evolutions arc obtained (Theorems 4.6, 4.9). The Bethe eigenvalue is shown to be 
a root of unity related to the generic period (Proposition 4.11). These results are 
derived and understood most naturally from the intrinsic picture on the dynamics 
as a straight motion in the 'ultradiscrete Jacobi variety' (4.30). We expect that the 
essential features explored in this paper persist in the generalized periodic box-ball 
systems [1, 2]. 

Acknowledgments This work is partially supported by Grand-in-Aid for Scien- 
tific Research JSPS No.l5540363. 

Appendix A. KKR bijection 

Put B = Bi^ ® ■■■ ® Bi^. The set V of our periodic box-ball system (2.7) 
corresponds to the B with the choice Zi = • • • = = 1. Let 

(A.l) V+ = {p&B\ew = Q} 
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be the set of highest elements. For the array of nonnegative integers m = {rrij) = 
(mi, m2, . . .), we put H = {j G Z>i | ruj > 0}. The data m is called a configuration 
if pj > for all j £ H, where 



(A.2) 



Pj 



L 

E 

i=l 



min(j, li)-2^ niin(j, k)mk 

k>l 



is called the vacancy number. The set H is necessarily finite and we parameterize 
it as if = {ji < ••• < js}- The data m is identified with the Young diagram 
containing rrij rows of length j, i.e., the rUj x j rectangular block for each j e H. 
We let Ad denote the set of all the configurations m. These definitions agree with 
the earlier ones for V+ (2.7), M (3.1), Pj (3.7) and H (3.8) when = ■ ■ ■ = 1^ = 1. 
Define 
(A.3) 



Rig(m) = {{Jl'^)l<^<^,,JeH G Z^^i X • • • X Z™- | < <■■■< < p,}. 

When Zi = • • • = Zl = 1 and the dependence on L is important, we write Rig^(m) 
in Appendices C and D. 

The purpose of this appendix is to describe the Kerov-Kirillov-Reshetikhin (KKR) 
bijection (j) = [5, 6]: 

(A.4) <i)-V+^ UmeM{{m, J)\J& Rig(m)}. 

An element of Rig(m) is called rigging. Each value jf"*^ is also said rigging. 
Similarly pj — is called co-rigging. The combined data (to, J) is called the 
rigged configuration. It is customary to depict it as the Young diagram m with 
each row attached with the rigging. We often exhibit the vacancy number pj on 
the left of the block of width j. The' riggings are arranged so as to decrease weakly 
downward within a block. Note that these definitions depend on the choice of B 
although we do not exhibit it explicitly. The choice other than li = ■ ■ ■ = li^ = 1 \s 
needed only in Proposition A.4. 

Example A.l. For the highest paths pi,P2 and ps in Example 3.2, one has B = 

Bf 19 and 



1 



4>{P2) = 



1 







The original KKR bijection [5, 6] is the one between the rigged configurations 
and the Littlewood-Richardson tableaux. The bijection (A.4) is obtained through 
a simple transformation of the Littlewood-Richardson tableaux and the highest 
paths [9]. Here we illustrate </) and casually along two examples rather than 
the systematic description which is already available in [5, 6, 28, 19]. Our convention 
here is opposite from [5] in the role of the rigging and co-rigging, and opposite from 
[19] in the order of tensor product. 

Regard a rigged configuration as a multi set of the pairs (row length, attached rigging). 
For example, the leftmost one in Example A.l is regarded as {(3, 1), (2, 1), (2, 0), (1, 8), (1, 4)}. 
Each element [j, a) of the rigged configuration is called a string with length j and 
rigging a. A string is singular if the co-rigging is zero, namely a = pj, which is the 
maximum allowed value in (A.3). We first illustrate the map (^~i. 
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Example A. 2. We begin with the basic case B = Bf , where the vacancy number 
(A. 2) reduces to (3.7). Consider the configuration m = (2, 1) for example. For 
L = 8, there are 6 rigged configurations depicted in the leftmost column of the 
following: 





21—1 


21—1 



10 2. 



10 2. 




21—1 



10 2. 



10 2^ 



1 

iro 



_|1 
1 



10 2^ imo 1^ omo 1^ ipi 

□o 2H0 ^do 



Ml A 

1 



i m o omo 1^ ipi 
3ni 2ni ^Ui 



^00 



_2^ IQO 0n0_2^ 




2 , 2^0 IDO 0n0_2, 



2, 2^0 IDO 0n0_2^ 



2, 201 IQI ^ 



_o ^ irno orno ^ imo i . omo 2 , iqi 2, a 1, a 

2 3n2 2Q2 

2 

The procedure to obtain the highest paths by applying (p^^ has been shown. Read- 
ing the numbers on the arrows backward, we find the image of those rigged config- 
urations under as follows: 



I(g)2(g)l(g)2(g)l0l02(g)2, 
I02(g)l(g)l(g)2(g)l(g)2(g)2, 

l(g)2(X'l(8'l(X'2(X'2(g)l(X'2, 
l(g)l(g)2(g)l(8)2(g)l(g)2(g)2, 
l(g)l(X'2(g)l(g)2(g)2(g)l(g)2, 
l(g)l(g)2(g)2(g)l(g)2(g)l(g)2. 



All these paths satisfy the highest condition (2.8). The KKR algorithm for obtain- 

"lL-i(rc')«)awithae Bi ={1,2}. 



mg 



^ proceeds recursively as 0^^^ (rc) 



This relation is depicted as rc rc' in the above. We have a = 1 G Bi and 
rc' = rc if the rigged configuration rc is free from singular strings. If there exist 
singular strings in rc, we set n = 2 G -Bi. In that case the new rigged configuration 
rc' is obtained by replacing any one of the shortest singular string (j, a = pj) by 
{j — l,Pj_i), where p'j-i is the vacancy number (3.7) with L replaced by L — 1 and 
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nii replaced rrij — 6ij + Sij-i. Namely one removes a box from the shortest singu- 
lar string and assings a new rigging to the shortened one so that it again becomes 
singular in the new environment. (If j — 1, just eliminate it.) Even when rc = rc', 
the vacancy number (3.7) is lowered by one by the change Bf^ Bf^-\ So one 
must revise pj in each step and keep track of L, which we did on the top line. 

To apply our inverse scattering method, the description in Example A. 2 suf- 
fices for In order to cover the content of Proposition A. 4 which is used in 
Proposition 3.4 and Lemma D.l, one needs to go beyond B = Bf^. We explain it 
along 

Example A.3. Take B = B2 ® Bi ® B2 ® ® Bi. We show the procedure for 
obtaining a highest path in B. 





1*1 




























0- 



iDo 



iDo 



000 ^ ^ H 



This time we have drawn a pair of diagrams in each step. The right ones are the 
rigged configurations. The left ones keep track of the 'shape' of the paths. For 
example in the top left diagram, the list of its row lengths 2, 1,2,3, 1 encodes the 
indices vn B = B2 ® Bi ® B2 ® B^ ® Bi where we start from. They are removed 
one by one from the bottom right as marked with •. The process corresponds to 
the reduction of the path shape as 

B2 (S) Bi (g) B2 (S) Bs (g) Bi B2 <8) Bi (g) B2 <E) B-i ^ B2 (S) Bi (g) B2 <8) B2 ^ 

B2 Si ^ B2 ^ Si ^ 0. 

In each step the vacancy numbers are revised by adopting the intermediate path 
shape as {li} in (A. 2). When a box (marked with •) is removed from the k th 
column, in accordance with (• • • ) S^ (• • • ) Bk-i, the rigged configuration is 

left unchanged if there is no singular string with length > fc. In this case we proceed 

1 .... 2 

to the next step by If there exist singular strings with length > fc, we have 

and remove a box (as marked with ★) from any one of the shortest such string. The 
shortened new string shall be attached with the rigging so that it becomes singular 
in the new environment. The algorithm ends up with the pair of 0. Reading 
the letters on the arrows backward, we get the sequence 11|2|11|221|2, where the 
symbol | separates those letters coming from different rows in the shape diagrams. 
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Reversing the letters within | • | , we find the image under <pg asll(8)2(g)ll(g)122(g)2e 
B, which is highest. 

The algorithm illustrated in Example A. 3 reduces to the simpler one in Example 
A. 2 when h = ■ ■ ■ = II = I- 

We have seen that the KKR algorithm for (j)^^ is a removal process of a rigged 
configuration creating a highest path from its rightmost component. Naturally, 
the map ({)B is an addition process building a rigged configuration by using the 
information of a highest path from its leftmost component. The rule of addition is 
easily inferred by looking at Example A. 3 backward. One regards the given path 
11 (g) 2 (g) 11 (g) 122 (Xi 2 € -B as the word 11|2|11|221|2 and adds a box (•) carrying 
these letters one by one from the left to form the prescribed shape diagram having 
the row lengths 2,1,2,3,2 from the top. If the letter is 1, one does nothing on 
the rigged configuration. If the letter is 2 and the change of the path shape is 
{■ ■ ■) Bk {■■■)(§> Bk+i, one adds a box (★) to any one of the longest singular 
string {j,a — pj) among those j > k. The new string of length j + 1 should be 
assigned with the maximal rigging so as to become singular in the new environment. 
If there is no such strings, one creates a singular string of length 1. 

For the highest path p g "P = Bf^ such that (j>'p{p) = ('n, J), it is easy to see 

n 

^^^^^^^^^^^^ 

<P-pig,Bf" (p 1 • • • 1) = {m, J). 

Proposition A. 4. Let p = hi® ■■■ &V = Bf^ be a highest path such that 
(j){p) = (to, J). Forn sufficiently large, define € V<S)Bf"' by u;(g)(p(g)l®") ~ ^(gu;. 
Then 4'-p0B'^"iO ~ {'^^^)> where the rigging I = {I'f ^) is given by /p-* = J^'^^ + 
min(/,j). 

In ^, sufficiently many components on the right arc also 1 G _Bi. 

Proof. We first show that (j)g^^j,^gf>,m{iii gi' p g:' 1®") = (m,I). In fact, the change 
from to4>g ^-p^B^" increases the vacancy number (A. 2) to pj+mm{l, f). 

Exactly the same increment has occurred from the rigging J to /. Therefore the 
co-rigging of (to, J) for V g) _B®" and the co-rigging of (to, J) for Bi (^V ® Sf " 
coincide, leading to 4'Bi®v®Bf'^ (w; g)pg) 1*^") = (to, /). On the other hand we have 

^V(»B'^"(»Bi ® by Lemma 8.5 in [28]. In view of 
the KKR algorithm, the last one is equal to ipp^Bf^iO- ^ 

The image 

(j>BHim,J)) ^ {h - xi,xi) (E) ■ ■ ■ (g) (II - xl,xl) E B ^ Bi^ (g) ■ ■ ■ (g) Bi^ 
can also be described by the picccwisc linear formula: 
ri(n) - Ti{n - 1) - ro(n) -|- ro(n - 1), 

n 

^^^^^(^ min(Zg, Uk) + {i - 2)z/fc - 4) - 2 ^ mm{uj , Vk)} , 

k g=l j<k 

where the maximum is taken over the subset v = {(vi, Ii) , (1/2, h) , ■ ■ ■} of the rigged 
configuration (m, J) regarded as the multi set of (row length, attached rigging). 
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(A.5) 



Ti(n) = 



Appendix B. Proof of Proposition 3.4 



Since Ei{uj{p)) ~ Ei{p) by Proposition 2.3, we may assume that wt(p) > 0. Such 
a path can be expressed as p = Tf{p+) for sonic d G Z and p+ G V due to Lemma 

3.1. Then Ei{p) = Ei{Tf{p^)) = Ei{p+), where the last equahty is due to Theorem 

2.2. Therefore it suffices to show (3.3) for p e P+, which we shall assume in the 
sequel. 

For any element of the form b & B = Bi^ ® • • • (g) Bi^ we define the quantity 
Di{b) e Z>o by 

C°wi 6 ~ 6* ® C^'(^)w; (b* e B, w e Bi) 

under the isomorphism Aff(i?i) ® Aff(i?) ~ Aff(B) ® A.E{Bi). (We have omitted 
the spectral parameters for b,b* .) When p is highest, its energy Ei{p) emerges not 
only in (2.10) but already in the relation (2.11) that produces vi. Namely, 

Lemma B.l. For any highest path p G V+, one has Di{p) = Ei{p). 

The following proof is direct but not intrinsic. 

Proof. Fixing a highest path = 6i (g) • • • (g) 6l, we regard -D;(p) and 2/2 in Figure 
2.2 as functions of there. The quantity Di{p) is the number of vertices of the 
bottom right type in Figure 2.1. Call them scoring vertices. Wc keep track of the 
scoring vertices appearing in Figure 2.2 as X2 is increased from to c = 2/2(0 ^ 
Di{p) = Ei{p) is shown if the scoring vertices for X2 = remain scoring and no 
new scoring vertices arc created during the increment of .T2. To sec this, recall the 
property 2/2 (2:2 + 1) = 2/2 (2^2) or 2/2(2:2 + 1) = 2/2(2^2) + 1 as noted in the proof of 
the Proposition 2.1. Since 61, . . . , 61, are fixed, a little inspection of Figure 2.1 tells 
that the increment of X2 never transforms the other type of vertices into scoring 
ones. On the other hand a scoring vertex may change into the top right type in 
Figure 2.1 if a = I — 1. Suppose this firstly happened as X2 = r is increased to 
X2 = r + 1 for some r < c(< I). Consider the leftmost such vertex that has ceased 
to be scoring at a;2 = r + 1. Such a situation is realized only if all the vertices 
(at a;2 = r + 1) on its left are bottom two types in Figure 2.1. Let a and /3 be 
the number of the bottom left types and the bottom right (scoring) types in them. 
Then one has r + 1 — a + P = I. On the other hand from the highest path condition 
(2.8), one also has a > /? + 1, where the last +1 is the contribution of the very 
vertex that has ceased to be scoring. Thus we obtain r = l + a — (3— 1>1, which 
is a contradiction. □ 

Proof of Proposition 3.4- Let p G V+ be a highest path such that (j){p) = (m, J) 

n 

hence /^(p) = m. Suppose that uf^ (g)p(g)l(g)---(g)l~p*(g)l(g)---(g)l(g) uf^ is 
valid. Hero wc take 1 <C A <C n. In general p* is a highest path longer than L 
(but much shorter than n). By Lemma B.l and H{ui (g Uk) = 0, we have Ei{p) = 

Di{p) = A(«?^opoio---oi) = A(p*oio---oiowf^) = A(p*oio---oi)- 

By Proposition A. 4, we know that (f){p* (g 1 (g • • • (g 1) = {m,I), where / = (/p^) 

reads /^'^ = J^'^ + A min(a, j). Therefore by taking A and a sufficiently large, one 

can achieve the situation !<•••< < • • • < /m] < /P"^^^ < • • • < < 
• • • . For such a rigged configuration {m,I), the KKR algorithm produces the path 
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p* (g) 1 (g) • • • (g) 1 = (j)-^{{m,I)) of the form: 



^3 

(B.l) . . . 2®-' . . . 2®-' . . . 2®-' s 2®-'+-'^ . . . 2®-'^"'^ . . . 2®-'^"^ . . . 2®-'^"^ 

Here . . . means an array of 1 € Bi and 2®^ G Bf-' . The number of appearance of 2®-' 
is nij. One can satisfy Lj ^ I for all Lj by taking A large. Moreover from the KKR 
algorithm, there are at least 1®-' between any two 2®-' . For such a pathp*0l0- • -01, 
it is straightforward to check Di{p* (g) 1 • • • 1) = m.m{l,j)mj. □ 

Remark B.2. In the pattern like (B.l), 2®-' is a soliton with length (amplitude) j. 
In this context the data m = {rrij) tells that there are mj solitons with length j. 



Appendix C. Proof of Proposition 3.7 

Lemma C.l. Let q e B®** be a highest path of length d and r e Bf^~'^ he a highest 

path of length L — d. Suppose that their rigged configurations are 4){q) — (1,1) and 
(j){r) = {n,K). Then the rigged, configuration of the highest path q ® r & Bf^ is 
given by (j){q (g r) = (Z U n, / U K'), where K' = [K'^ "'^) is given by 

(C.l) K[^'^ = K^^ - 2 ^ min(j, k)k + d, 

k 

where {lUn,I U K') means the union regarding {1,1) and {n,K') as multi-sets of 
rows assigned with rigging. 

Proof. The new vacancy number p'^ for the width j block in {lU n,I U K') reads 
p'j = L — 2 2^ min(j, k)lk — 2 min(j, k)nk. When applying the KKR map (t)~^ 
to (Z U n, 7 U K'), the co-rigging of the i-th row of the width j block in (n, K') is 

(C.2) p'^ - K[^'^ = iL-d-2j2 niin(j, k)nk) - K^^l 

k 

The right hand side is equal to the co-rigging of the same row in the rigged config- 
uration (n,K). Meanwhile the co-rigging of the {1,1) part is not less than that in 
the original q plus the vacancy number of r. Therefore the algorithm of the map 
(j)-'^ proceeds as (^"^((Z U n, 7 U if')) = <;^"^((Z,7)) r = g r. □ 

Example C.2. Take q = 1112122, r = 111221221122, hence d = 7,L = 19. The 
highest path p^ in Examples 3.2 and A.l is expressed as p2 = q®r. 



m = 



(C.3) 







|o 






2 




3 





The rigged configuration for 4'{p2) in Example A.l contains this (/)(g) indeed. The 
rest of it corresponding to {n, K') is related to c6(r) = (n, K) as 







□ 1 






|0-2 X 3-1-7 






3 






2-2x3-^7 




6 






3 


-2x2+7 



in agreement with (C.l). 
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Lemma C.3. Let q G .Bf"' be a highest path of length d and r G Bf^~'^ he a highest 
path of length L — d. Set (j){r (g) g) = (m, J) and 4>{q (g) r) = (m', J'). Then m = m' 
and i{J') — i{J) + d G J{m) are valid. 

See (3.12) and (3.16) for the definitions of ~ and l. 

Proof. We use the same notation as in Lemma C.l. Thus m' = ZUn and J' = lUK'. 
Applying Lemma C.l with {I, I) and (n, K) interchanged, we find m = l[Jn = m' 
and J = I'yjK with = — 2 niin(j, k)nk + L — d. In terms of the vacancy 
number p'^ for m in the proof of Lemma C.l, the relation between the rigging J 
and J' is summarized as 

(C.4) = + 2 ^ min(j, k}lk + p', - d, 

k 

(C.5) iiTp) = K.^'^ + 2 ^ min(i, A;)/fe - d. 

k 

See the picture after (C.8). Now switch to the extended sequences l{J) and i{J') 
by (3.16). Apart from — d, this is exactly the effect of the shde Ofc ^1'° ("^-H) 
the width j block of i{J'). Therefore we conclude b{J) = {JJ^ cF^)b{J') — d. □ 

Proof of in Proposition 3. 7. Without loss of generality we may assume d' = 
and < d < L. Since p+ and ^'^^ highest paths, they must have the form 
p+ = r ® q and p'_^_ = q®r, where the both paths q e Bf'^ and r e Bf^~'^ are 
highest ones. Then the assertion follows from Lemma C.3. □ 

To show in Proposition 3.7, we explore the full implication of the right hand 
side of (3.17). Let H = {ji < ■ ■■ < js} be the set of lengths of rows of m as in 
(3.8). 

Lemma C.4. Suppose l{J) + e ~ i.(J') + e' G Jim) for some m = {mj) G M., 
J, J' G Rig(m) and e, e' G Z. Then there exist < Ij < nij {j G H) and < d < L 

such that the relation (tI' ■■■(Tj;t{J') = t{J) + d holds. 

Proof. Prom the assumption, there exist Ij's and d such that aj^^ ■ ■ ■ a-^ b{J') = 
l{.J) + d. Our task is to show that they can always be chosen within the range 
< Ij < nij {j G H) and < d < L. From the last comment in Section 3.2, to 

achieve the relation a^^^ ■ ■ ■ a^" i(J') = l{J) + d by the successive transformations 

t(J') ^ a}-i(J') ^ 'T^lr>j7t(J') ^ > • • •fTj-t(J'), one must apply a^' 

so that the blocks i/j^ , J7jfc+i , • • • , JTj, in CTj^^.'' ' " '-i^') already coincides with 
that in l{J) up to an overall additive constant. Having this in mind we proceed to 
showing < Ij < rrij. 

Certainly, Ij^ can be taken as < Ij^ < ruj^ by using (3.13). Since the slides 
are invertible, one may also assume that the first non-zero number in the sequence 
lj^,ljs-i, • • • is positive, namely, Ij^ = Z,,_i = . . . = Ij^^^ = and Ij^ > 0. If there 
is no such t, it follows that l{J') = i{J) + d and we are done because d > may 
be assumed without loss of generality and then d < L is obvious from L > pj (3.7) 
and (A. 3). Henceforth we assume that Ij^ > exists for some 1 < t < s. We first 
claim that Ij^ < nij^. In fact, the case t = s is within our assumption. If Ij^ > ruj^ 
for t < s, the upper blocks jT/j (/3 > jt) acquire the uniform shift 2jtlj^ under a^* . 
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On the other hand the rigging j'^^^ gets shifted at least by 2pj^ + 2jtlj^ . Therefore 

to adjust <T-H{J') to some i{J), one must extract at least pj^ + 2jtlj^ from the 
new rigging on the upper blocks J'fj as an overall constant. But this fails since the 
original rigging there is not greater than pp{< PjJ hence the result of the above 
extraction is not greater than {pjs + 2jtlj^ ) — (pj^ + 2jilj^ ) < 0. Thus we have verified 

Next we prove < Ij^ < mj^ (1 < a < t) by induction on a assuming that 
< Ij, < ruj, for all a + 1 < 6 < t. Setting K = (kI^^) = a^^ ■ ■ ■ i{J'), we have 
i^l^'^ = J-^iil+'2j:Ua^^ijf,3c)lj.- Note that 



t 

(C.7) k[^'^ < Pj, +2j2^i^{js,jb)ln, 

b—a 

where the former follows from Ij^ > and the latter does from /j_ < m,^ and Jrn < 
Pj^. As explained above, the blocks J'jj (/3 > ja) in K — d should already coincide 
with those in l{J) for some d. To show Ij^ > 0, suppose Ij^ < on the contrary. 
Then we have xj^"^ < - pj^ + 2ELamin(ia,.ifc)Zj, < 2 ^Lo 
Prom this and (C.6), K — d can coincide with some l{J) only if 

t t 
2 X] ' ^b)ln <d<2^ min(ja , ■ 

b=a b=a 

But this is impossible because of Ij^ > and the induction assumption, verifying 
Ija ^ 0- To show Ij^ < rrij^, suppose Ij^ > nij^ on the contrary. Then we have 

Kmjl > 2pj^ + '^J2b=a^^^ija' jb)ljb- ^om. this and (C.7), -fC — Jean coincide with 
some t{J) only if 



t t 



^'ia + 2 X! ^^^ijo,,jb)lh < J < Pi» + 2 E mm{j„jb)lj 

b=a b=a 

Again it is easy to check that this is impossible by the reason similar to the previous 

inequality, proving Ij^ < ruj^ . 

Finally we consider d in the relation L(J)+d =- a^^ ■ ■ ■ o"-^ i^{J') — <Jj^ ■ ■ ■ <Jj' t-{J')- 
Prom (C.6) and (C.7) with a = 1, we get Jm^ + d > pj^ +2 J2b=i ^^^Ut, jb)ljb and 
+ d< pj, + 2X;Li ^Mjs,jb)lh- Therefore in order that (jf^^) e Rig(m) to 
hold, d must satisfy < 2J2l^^ mm{jt,jb)lj^ < d < pj^ + 2J2l^i min(is, < 
Pje + 2X^5=1 min(js,j(,)TOj^ = L. In particular, d = can happen only if VZj^ = 
0. □ 

Proof of in Proposition 3. 7. By Lemma C.4, J and J' are connected by the 
relation (life <''fc')''("^') = ^-iJ) + <^ some < < rrij and < rf < L. If 
VZfe = 0, the assertion follows easily from the definition of the KKR bijection 4>. 
Henceforth we assume that Ij > for some j G H. Let n = {rij) and I = (Ij ) be 
the Young diagrams giving the decomposition m = I L) n, where < rij < mj is 
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defined by rij = nij — Ij. By the last remark in tlie proof of Lemma C.4 we have 
< < L. 

Define /(^^ = (4''^)i<i<(, and i^(^) = {Kl'^)i<i<nj by 



(C.8) (/}• 



0) 



'0) 



Then define further il'i^^^)i<i<ij and iKy')i<i<nj by (C.4) and (C.5), respectively. 

By construction we have J = I' U K and J' = I U K', where {Jn-'+i^ ■ ■ ■ ■> Jm]) = 

. . . , and ( fffi, • . • , J^^) = {K'l'\- • • , Klip). The following figure is 
helpful to grasp these relations. 



t 



K'(j) 



jU) 



t 



' t 



T'U) 



jU) 



We claim that I e Rigd(/) and K e Rigi_rf(n). In fact, < /p^ < ■ ■ < ll 



IS 



obvious and /, 



0) 



T'U) 



2 Y,k k)lk - p'j + d < d - 2Y,k k)lk which 

is the vacancy number for Bf'''. (Here p'^ = L — 2 min(j, A:)TOfe is the same as 
that in the proof of Lemmas C.l and C.3.) Similarly, < k[-'^ < • • < Kn] is 



obvious and = K^^' + 2 min(j, A;)^^ - d < p'^ + 2 min(j, A;)//c - d = 

L — d — 2 min(j, k)nk which is the vacancy number for Bf^~'^. 

Set q = e and r = <?!)-i((n, X)) € Sf^"''. Then Lemma C.l 

tells that the highest paths p+ and p'_^ in Proposition 3.7 are given by = 
(j)~^{{in, J)) = r ® q and = 4'~^{{^^J')) = q ® r. Therefore we obtain 
= Tt{p+). □ 



The proof of Proposition 3.7 is finished. 



Appendix D. Proof of Theorem 3.12 

For distinction we write the time evolution of the angle variable (3.14) as r; 
within this appendix. Obviously TiTk = TkTi is valid. 

It is straightforward to check the commutativity of the diagram (3.18) for I = 1. 
In fact, if p = Tf{p+) for p+ G V+{m) with 4'{p+) = {rn,J), we have $(Ti(p)) = 
<^{T^+^{p+)) = [l{J) +d+l]. On the other hand due to ri(i(J)) = l{J) + 1, we 
obtain ri($(p)) = n{[i{J) + d]) = [l{J) + d+l]. 

The commutativity ti$ = ^Ti and Lemma 3.6 reduce the proof of Theorem 3.12 
to the highest paths p+ G V+{m). In fact, the equality $Ti(p) = Tj$(p) for general 
pathp = T|*(p+) e ■P(m) is deduced from $T;(p+) = Ti^{p+) by multiplying rf on 
the both sides and using the commutativity TiT; = TiT/ and tiT; = tjTi. 

In the remainder of this appendix we fix the Young diagram m = {nij) G ^A and 
assume that p G 'P+(m) C Bf^. {rrij is the number of the length j rows in m as in 
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the main text.) Then we determine J, K, I, d and q successively as follows: 
(D.l) (t>{p) = {m,J), JGRigi(m), 

(D.2) Tii{J)=K &J{m), 

(D.3) [K] = [(.(/) + d\€ J{m), I e Rigi(m), d>0, 

(D.4) ct>-\{m,I))=qGV+{m), 

where the definition of t is available in (3.16). The role of these objects will be 
seen clearly in (D.8). In (D.3), the choice of / and d is not unique. However what 
matters in the following proof is the combination Tf{q) whose uniqueness has been 
assured by Proposition 3.7. 

For N denote by — (Nrrij) the iV-fold repetition of m and define the 

map 

(D.5) 

TTAT : J(m) — ^ Z^™^i X • • • X Z^™^» 

({.J\'%^-L, . . . , (#^)iez) - ((jF^^)i<i<iv™,, , . . . , {J\''\<i<^^i}j ■ 
In what follows we set B = Bf^ and assume that C > {N + 1)L. 
Lemma D.l. For 'jtn{l{J)) G Rig j^i^{m'^), one has 

N 



(D.6) 0B^((m^,7rjv(t(J)))) =pOpOp(S) OpOl O 1 

Af-l 



(D.7) (/.^'((m^, 7rAr(X))) = p* (g) Ti{p) (g) Tiip) ®---® Ti{p) ®p** ® I ® ■ ■ ■ ® I, 
where p* € Bf^ is defined by (2.11) and vi ® l®'^ ~ p** ® ui. 

Proof. The relation (D.6) is derived by repeated use of Lemma C.l. To see (D.7), 
note from (2.11) that 

N N-1 



0p0 • • -0^01 • • • O 1 (g)T;(p) • • • (8 Tlip) Op** 1 • • • O 1 O U;. 

From this and (D.6), the relation (D.7) is obtained by applying Proposition A. 4 
with p and (m, J) replaced with p'®^ and (m^, 7rAr(/,( J))). □ 

Take the slide cr = ■ ■ ■ a^^ e A that achieves (D.3) via (7{K) = t(/) + d. By 
the argument similar to the proof of Lemma C.4 (especially after (C.7)), one can 
show that Ij's can be chosen as < Ij < rrij. We introduce the Young diagram 
rh^ = {rhf) by = Nrnj — Ij and let ttjv be the map (D.5) with Nrnj replaced 
by m^. 

Lemma D.2. As C andN grow large satisfying C > NL :$> 1, thepath(j>'^^{{fh^ ,TTN{a{K)))) 
takes the form ^ ?7, where the left part ^ is independent of C, N and the right part 
r] growing with C,N is the same as the corresponding part in (l)~^{{m^ ,-kn{{K))) 
(D.7). 

Proof. From the definition (3.11) and a = a-)- • • • cr^^^ , the rigged configuration 
(fh^ ,ttn{o'{K))) is obtained from {m^ ,nN{K)) by removing the bottom Ij rows 
and adding 2 J2b=i ^^^U' jb)hb to the rigging for each block of length j rows. This 
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operation does not change the co-rigging of the remaining part. Therefore the 
assertion follows from, the KKR algorithm. □ 



Proof of Theorem 3.12. We are to check = in the following diagram: 



P I- 



n 



(D.8) 



Tiip)=Tf{q) 



H [K] = [c{I)+d] 



Consider the rigged configuration {rh^ ,TtN{cr{K))) = (m^, 7TAr(t(/) + d)). The 
latter expression tells that this is also obtained from (m^, 7rAr(t(/) + d)) by removing 
the top Ij rows in each block of length j rows. Thus from the KKR algorithm, the 
paths (f)'^^{{m^ ,TtN{i.{I) + djj) and (f)'^{{rn'^ + d))) are the same except 
only some components that are located on the right and remain unchanged with 
growing C, N . On the other hand from the definition (D.4) and Lemma D.l (D.6) 
with (p, J) replaced by (g, /), we find 



N 



</>Bl((TO^,^ArW/)+(i))) = l^ 



1. 



Combining the foregoing argument and Lemma D.2, we can equate the middle part 
of this with that in (D.7), leading to T;(p) = Tf{q). □ 
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